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, Two-dimensional interacting electron systems become strongly correlated if the electrons are sub- 

. ject to a perpendicular high magnetic field. After introducing the physics of the quantum Hall 

regime the incompressible many-particle ground state and its excitations are studied in detail at 
' fractional filling factors for spin-polarized electrons. The spin degree of freedom whose importance 

was shown in recent experiments is considered by studying the thermodynamics at filling factor one 
and near one. 
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I. INTRODUCTION 
A. Overview 



Two-dimensional electron systems (2DES) of high-mobility can be experimentally realized in modern semiconductor 
*~1 ■ devices near interfaces. In particular, the limit of a strong magnetic field normal to the interface offers the opportunity 
"j^ ' to study an electronic system, which is unique in many respects. Due to the strong field, all electrons can be 
: accommodated in the lowest orbital Landau level of quenched kinetic energy resulting in a macroscopic degeneracy of 
'^l^ the non-interacting ground state. Thus, the only relevant energy scale of the system is set by the interaction energy, 
'■^ which turns the system into a strongly correlated one. 
^ The existence of energy gaps at certain filling factors is one of the reasons for remarkable anomalies in magneto- 

O I transport, the quantum Hall effects. These exhibit features that are very different from those in lower magnetic field 
experiments. In the case of small occupancy of the lowest Landau level, the existence of gaps can be traced back to 
the occurrence of an incompressible many-particle ground state. 

In this Introduction, we review the most important experimental and theoretical developments of the physics in the 
^ ■ quantum Hall regime and derive the appropriate many-particle Hamiltonian. Then, we discuss the problem of two 
interacting particles of equal and opposite charge, which frequently appears in the many-particle treatment. 

In the second part, we study in detail the properties of the many-particle ground state and its elementary charged 
' excitations, the quasiparticles, in the limit of an infinitely strong magnetic field neglecting the spin degree of freedom. 
This is explicitly done on the basis of a special model and by means of numerical methods in the case that only one 
PsJ _ third of the available one-particle states in the lowest Landau level is occupied, i. e. at filling factor 1/3. Furthermore, 
. we discuss the generalizations to filling factors of the form p/q, where p, q are integers and q odd. 

The third part is devoted to the consideration of the spin degree of freedom, which recently received a lot of attention 
from experiment. It is shown that at certain filling factors the ground state exhibits spontaneous spin magnetization, 
y which defines the system as a quantum Hall ferromagnet. After the discussion of the ground state properties near 
filling one, where novel charged spin texture excitations appear, a many-particle theory is developed to determine 
(— ( thermodynamic properties at exactly filling factor v = I. The theory accounts for the spin- wave like excitations above 
Q ' a completely spin-polarized ground state. The many-particle theory improving the inadequate Hartree-Fock theory 
O ■ predicts the temperature dependence of the spin magnetization. The results are compared with recent experimental 
^ ' data as well as with different theoretical approaches. 
• ^ . Our considerations are based on a microscopic electronic Hamiltonian. It becomes clear that this strongly correlated 
' 2DES shares many difficulties with theories of metallic fcrromagnetism, yet it is free from the consequences of a complex 
^ , bandstructure, which have confounded comparison with experiment. 
■ - - I Eventually, we speculate about possible routes of future research in this lively field that has offered surprises for 
more than twenty years. 
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B. Two-dimensional electron systems in a magnetic field 



Traditionally, any microscopic description in solid state physics starts from simplified models trying to capture the 
essential features before systematically extending the model in order to reconcile it with the often much more complex 
experimental results. Pursuing such an approach, we describe electrons of spin S = 1/2 and charge q = —e < 
moving in a three-dimensional crystal, which are subject to a spatially constant and time independent magnetic field 
B = Be^ = B±ez (-B > 0), in a first approximation as independent particles in a continuum. Due to the separation 
of the motion in the x-y-plane and the z-direction, the Hamiltonian can be written in Cartesian coordinates as 

Ho = ^Jp-qAr + ^^~g.,,BS. (1) 

Here, the kinematical momentum in the plane is given by p = {(h/i)dx — qA^, {h/i)dy~qAy), whereas the momentum 
Pz — {^/'i)dz describes the free motion in the z-direction. The gauge freedom of the vector potential allows to write 
it as A(r) = B{—ay, (1 — a)x,0) with the parameter a G R. The gauge parameter a covers for a = and a = 1/2 
the two mostly used gauges, viz. the Landau gauge and the symmetric gauge. The last term is the Zeeman energy 
and is due to the coupling of the dimensionless electron's spin S to the magnetic field {ge - gyromagnetic factor, /i^ 
- Bohr magneton of the electron) . The natural length scale is set by the cyclotron or magnetic length £c — y/^7\qB\, 
where ic[nm] = 25.65/ y/B[T], which exceeds for accessible fields up to i? ~ 20 T the lattice constant by more than 
a factor of 10 and justifies the application of a continuum model. 

The solution of Eq. (||) separates in three simple energy eigenvalue problems: first, the motion in z-direction for 
finite extension w with energies Ez^i = j/(2rne), where kz^i oc 1/w as long as there is no in-plane component of 
the magnetic field, second, the spin degree can be simply described by the two eigenvalues of Sz = ±1/2, which are 
denoted by the majority spin direction "f for a = +1 and the minority spin direction | for cr = — 1 as long there is 
no spin-orbit coupling, and third, the solution of electrons in the plane that can be reduced to a harmonic oscillator 
problem, whose energy scale is hujc, and where ujc = \qB\/mf, is the cyclotron frequency. 

The problem becomes effectively two-dimensional if one is able to create well distinguished energy levels, so-called 
subbands, by requiring that the temperature T is smaller than the energy difference Aii^oi between the first excited 
and the lowest subband energy, i. e. the following inequality holds 

keT < AEo, = . 2 

The estimate is based on the assumption of a rectangular quantum well of infinite height in z-direction. In order to 
satisfy such a condition for, say, temperatures below 200 K, a width smaller than at least 120 A is necessary. Only 
the advent of modern semiconductor technology made the practical realization of devices possible, where potential 
changes within such a short distance could be achieved. By means of appropriate doping, it is possible to create a 
two-dimensional electron system (2DES) at the semiconductor-oxide interface whose Fermi energy can be controlled 
by a gate voltage. Since the end of the seventies, technological progress made the usage of a new class of semiconductor 
devices for such experiments possible, viz. GaAs-Al^Gai-xAs heterostructures. They are based on the appropriate 
doping of stacked GaAs and Al^Gai-xAs layers, which form either heterojunctions in case of one interface or one 
and many quantum wells, respectively. In particular, the growth of heterostructures by means of molecular beam 
epitaxy allows a much more controlled manufacturing, which leads to a reduced number of impurities and, hence, to 
an increasing mobility of the samples. 

Assuming that the electrons occupy at sufficiently small temperatures only the lowest electronic subband, the 
one-particle problem neglecting the subband energy contribution has degenerate eigenvalues 

En,a = ?it^c (^n + ^ ' (3) 

which are labeled by the orbital quantum number n — 0, 1,2, . . . denoting the Landau levels and the eigenvalue of 
the z-component of the spin operator, while huJc determines the distance between two adjacent Landau levels [0. 
Note that the situation becomes different in a real 2DES, since one has to account for the effective values of the 
mass, gyromagnetic factor, and dielectric constant of the host semiconductor. For GaAs, where the effective mass is 
771 — 0.067me, the gyromagnetic factor g — —0.44 [^, and the dielectric constant e = 12.7eo, the energy spectrum is 
schematically depicted in Fig. |l|. The explicit form of the eigenfunctions, belonging to the eigenvalues depends 
on the gauge. In symmetric gauge, the energy eigenvalue problem is reduced to the solution of a two-dimensional 
isotropic harmonic oscillator problem, and in Landau gauge, it describes a free particle subject to periodic boundary 
conditions in the y-direction and a one-dimensional harmonic oscillator in the x-direction. In either case, each of the 
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FIG. 1; The two lowest spin split orbital Landau levels as in a GaAs heterostructure, where < ^ huic- The degenerate 
one-particle states are denoted by circles. The ground state of non-interacting electrons with majority spin direction t for — 1 
is indicated by filled circles. 



Landau levels is degenerate with respect to a quantum number determined by the respective gauge. For the symmetric 
gauge, these are the eigenvalues of the z-component of the angular momentum = {h/i)d^, which can be labeled for 
each Landau level by m = — n, —n + 1, . . ., and in the Landau gauge by the momenta in the y-direction ky^i = 2^1/ Ly 
[Ly - width of a rectangle in the y-direction, I = 0, ±1, ±2, . . .). To be definite, let us quote the eigenfunctions of the 
lowest orbital Landau level (LLL) n = with m = 0, 1, 2, ... as a function of the dimensionless complex coordinate 
z = X — iy = {r^ — iry)/ic (note sign{qB) = —1) 

1 1^1^ 

<Pn=0.7n{z) = ^Pm{z) = . Z^e « . (4) 

v27r2™TO! 

The orbital degeneracy iV$ = $/$o of each Landau level counts how many elementary flux quanta $0 = h/e are 
contained in the total flux <i> = BA through the area A of the plane, i. e. Ng, = A/(27r£^). This implies for a 
system of finite extension in the plane that the degeneracy of a Landau level is finite too, e. g. on a circle of radius 
R, the degeneracy is i?^/(2£^). In case of the LLL, the eigenfunction (pm{z) oc covers an area of 27r(m -I- and 
encloses therefore (m -I- 1) elementary flux quanta. In other words, the power of z, which is also the number of zeros 
of the eigenfunction, equals approximately the number of flux quanta covering the area of the system. An important 
quantity is the filling factor v, which is defined as the ratio of the particle number N to the orbital degeneracy 7V$, 
i. e. 

V = iV/7V$ = 2Tilln , (5) 



where n = N jA is the particle density.] 199 1 We note that the degeneracy is proportional to the strength of the 
magnetic field B. This allows to accommodate an increasing number of electrons in each orbital Landau levels when 
the magnetic field becomes stronger, or for a fixed particle number, the Fermi energy is swept downwards with 
increasing magnetic field. If the field is so strong that at zero temperature all electrons reside in the lowest orbital 
Landau level, we will loosely speak about the strong field limit . 

Besides the symmetries, which are related to a specific gauge, and which are satisfied even for finite systems, there 
exists a continuous symmetry of the Hamiltonian Eq. (|l|) for an infinitely extended system since the generator t of the 
magnetic translation operator commutes with the Hamiltonian Hq, i. e. [t,_ffo] = 0. The generator of the magnetic 
translation in the plane is the canonical momentum t = p + gA, and the magnetic translation operator is hence given 

by i 

Tro = e-'^''°(P+'?^) . (6) 

Before closing, let us emphasize that such a model of independent electrons ignores at least two important ingre- 
dients, viz. disorder and interaction between the electrons, which become important when discussing the strong field 
limit. 



C. A short historical account of the 2DES in a strong magnetic field 



In this Subsection we review some of the remarkable experiments of a 2DES in a strong magnetic field, which have 
made the physics of strong magnetic fields such a lively field of research. Up to the discovery of the integer quantum 
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Hall effect (IQHE) , the physics of a 2DES did not seem particularly spectacular despite the new opportunities due to 
the reduced dimensionality. For example, the system behaves like an ordinary Fermi liquid without magnetic field 

The measurement of transport properties belongs to those methods that allow the phenomenological characterization 
of solids. Hall effect measurements serve as the experimental standard method to determine the sign and concentration 
of the charge carriers, and they are explained in the framework of the classical Hall effect. It is evident that already 
at the classical level an independent electron model without disorder is insufficient to describe transport. Sources 
of disorder in a semiconductor 2DES are charged impurities, e. g., donors, near the 2DES or the roughness of the 
insulator-semiconductor interface in metal oxide semiconductor field effect transistors (MOSFETs). In the Drude 
theory of the Hall effect, carriers of charge q move in the electric field E = E^Gx and cause a current density j = ctE 
((7 - conductivity tensor), which can be determined by inversion of the resistivity tensor tr = . Scattering of the 
charge carriers on impurities leads to elastic scattering events with elastic scattering time r, and the Drude theory 
gives for the components of the resistivity tensor p 



Po 



nq'^T 



Pxy 



h 
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(7) 



B _ 1 
nq V q 

where po is the resistivity of the Drude theory without magnetic field. Such a model is often sufficient to explain 
experiments in fields up to 0.1 T. 

An unprecedented development started with the observation of the IQHE in 1980 by von Klitzing when he studied 
the Hall resistance Rh = \pxy\ and the longitudinal resistance R oc pxx of a silicon MOSFET with a, at that time, 
comparably high mobility at temperature 1.5 K in high magnetic fields The Hall resistance exhibited plateau 
values Rh — {hle^)/v around integer filling factors = 1,2, . . ., while the longitudinal resistance started to vanish 
in the middle of the plateaus showing deep dips, c/. Fig. |2[ These features become most pronounced at very low 
temperatures, where the Hall resistance exhibits a staircase like structure. Simultaneously, the longitudinal resistance 
vanishes except for peak- like structures around v = (n -I- 1/2), where the transition from one Hall plateau to the next 
one occurs. Particularly remarkable is the accuracy of the Hall resistance in the plateau region, whose values are 
integer fractions of the von Klitzing constant Rk — hje^ . Currently, the relative standard uncertainty has a value of 
3.7 10~^ 0, and the IQHE serves as the metrological standard for the resistance. These results are clearly different 
from those of the Drude theory, since the Drude result for the Hall resistance merely agrees with the IQHE result at 
those where the linear p^y vs. B curve (0) crosses the staircase structure. This happens at those values for the 
magnetic field, where v — ft,n/(|ei3|) attains an integer value, c/. Fig. |^. 

When in 1982 Tsui et al. performed the same experiment with a Ga^s-sample of higher mobility at lower temper- 
atures in stronger fields, they found similar features in the Hall and longitudinal resistance around certain fractional 
filling factors of the form J/ = p/g with odd (7, e. g.,fovv = 1/3, 2/3, 2/5 to name the most prominent ones, see Fig. |^. 
Hence, the name fractional quantum Hall effect (FQHE) was coined M. From a phenomenological viewpoint, these 
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FIG. 2: Results of transport measurements in a high-mobility GaAs-heterostructure (/i ~ 1.3 x W^cm^ /Vs) at low temperature. 
The Hall resistance Rh = Pxy exhibits plateaus at integer (IQHE) and fractional (FQHE) filling factors u, while the longitudinal 
resistance R = pxxiLx/ Ly) shows deep dips in the plateau region indicating the vanishing of pxx when approaching zero 
temperature. From 



quantum Hall effects are quite similar, although the energy scales involved significantly differ. For the understanding 
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of the quantum Hall effects (QHEs) it is important to note that at those where the effects appear, the disorder-free 
system exhibits an energy gap above the ground state. In the FQHE, this is a gap above the many-particle ground 
state caused by the dominating electron-electron interaction, cf. Sect. ||, whereas in the IQHE, the gap has its origin 
in the energy gap between adjacent Landau levels at integer filling. Since the quantum-mechanical treatment of a 
disorder free system shows a linear B dependence as the Drude theory, a certain amount of disorder is necessary in 
order to create plateaus. On the other hand, sufhciently large disorder lets the QHEs disappear, at first, the FQHE, 
at larger disorder, the IQHE too. This can be explained by the different magnitudes of the energy gaps: the energy 
gap in the IQHE is approximately one order of magnitude greater than that in the FQHE. Hence higher mobility 
samples than for the IQHE are needed to observe the FQHE as larger disorder washes the energy gap out. A similar 
effect can be observed by increasing the temperature, which causes the disappearance of the FQHE, and at higher 
temperatures that of the IQHE |lOj . 

The theoretical explanation of the IQHE is currently based on a one-particle model with quenched disorder. In 
a two-dimensional system in a strong magnetic field, disorder causes the broadening of the Landau levels forming a 
continuum of states ||ll| . Most of the states are localized and only those states exactly at the energy of the originally 
discrete Landau level (LL) energies seem to be extended jl^. Such a system is an extensively studied example of a 
localization-delocalization transition in two dimensions. The formation of a mobility gap due to localized states has 
a similar impact on transport as a charge gap in a pure system [p^ . There are various ways to make this statement 
more precise. For example, using a Kubo formula for the Hall conductivity axy, one can show that only extended 
states below the Fermi energy contribute Hence, increasing the chemical potential up to the center of the 

LLL /i — huJc/2 keeps axy = 0. When reaching the half-filled LLL, the conductivity jumps to the value e^/h and 
remains again constant up to a filling with a Fermi energy at the center of the next Landau level. Other explanations 
are based on a gauge argument and the role of the edge of the system [|l5|, However, we are still far from a 
microscopic transport theory of the IQHE JTtI that makes quantitative predictions about the influence of disorder 
and temperature on the conductivities. Nonetheless, there is the belief that this explanation meets the essence of the 
effect. In the FQHE, the energy gap is not the result of a one-particle effect but is caused by electron correlation. 
Its origin can be traced back to the specific properties of the many-particle ground state at the so-called magic filling 
factors p/q. In particular, Laughlin's proposal for a ground state wavefunction at v = 1/q paved the way for a 
deeper understanding ||l^. One of our main subjects is to explain the existence of these gaps. The importance of 
electron-electron interaction is supported by the experimental fact that the FQHE features become more pronounced 
when increasing B, but keeping v constant, and hence diminishing the influence of higher LLs. We see that, contrary 
to the IQHE, even the study of the ground state and the low lying excitations of interacting electrons in strong fields 
is a difficult subject on its own. This is also one of the reasons why a lot of work was devoted to this issue, whereas 
the transport properties of the FQHE are usually explained by adopting arguments from the IQHE to the case of the 
FQHE 

Further inspection of Fig. ^ shows that pxx behaves around v — 1/2 quite differently. This is attributed to a 
Fermi liquid like ground state of weakly interacting particles, so-called composite fermions, at = 1/2. In particular, 
surface acoustic wave measurements gave first evidence of this effectively one-particle picture ^o| , Moreover, 
theoretical work p2| based on a singular gauge transformation supported such a viewpoint. 

A third group of experiments covers the very low filling factor region, where one expects a behavior typical for a 
Wigner crystal pinned by disorder p^ . Various experiments in the low filling factor region indicate that physical 
properties change when compared with those in the FQHE region. Transport experiments even at ~ 1/5 ||24| , e. g., 
show a diverging longitudinal resistivity indicating an insulating ground state and reentrant behavior in the typical 
FQHE region. So far, however, an experiment giving unequivocally evidence for the existence of such a ground state 
is still missing [ p3[ . 

Early experiments in the QHE region explored almost exclusively dynamical properties examining the response 
of the 2DES to electromagnetic fields. This was mostly done by studying transport properties. The investigation 
of the thermodynamic behavior was particularly pushed by predictions on novel spin structures, although Halperin 
had already noted very early the importance of the spin degree of freedom in Ga^s- heterostructures tea . The spin 
magnetization can be investigated, for example, by nuclear magnetic resonance (NMR) experiments pq], magneto- 



absorption experiments [pTj , and photo- luminescence measurements 28 



Transport experiments at higher LLs emphasized the importance of the spin, and &i v — 5/2, the only FQHE with 
even denominator in a one-layer system was found ^ . Recently, there were also unexpected observations in the higher 
Landau level range. Near half filling of higher LLs, e. g., at v = 9/2, anisotropies in pxx were found. They are caused 
by ground states with a strip like density distribution, which is different from the LLL behavior p^ . 

Within the last ten years, further experimental information about the physics of the quantum Hall regime was 



gathered by studying double layer systems, where two parallel quantum wells are brought close together |30, pll. In 



such systems, the phenomena become much richer due to the loss of symmetry of the Coulomb interaction and the 
possible tunneling of charge carriers from one layer to the other one. 
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Hallmarks of the theoretical understanding of the QH-regime are Laughlin's seminal proposal of a ground state 
wavefunction at z/ = 1/q, the introduction of the quasiparticle concept and the construction of a hierarchical theory 
to explain the occurrence of the FQHE at v = p/qhy Halperin and Haldane and Jain's composite fermion approach. 
The construction of excited states that describe gapped density like intra Landau level excitations above the Laughlin 
state allowed the prediction of a roton-like minimum . 

The construction of a field-theoretical approach was initiated by the deep observation of Girvin and MacDonald 
that the application of a special singular gauge transformation to the Laughlin state at 1/q by attaching q flux 
quanta per electron leads to a bosonic state ||33|| . Most remarkable, the resulting bosonic wavefunction exhibits in 
contrast to the original fermionic Laughlin state an algebraic off-diagonal long-range order (ODLRO) of the one- 
particle density matrix. This observation fostered attempts to search for an order parameter describing the unique 
Laughlin state and the derivation of a phcnomcnological Landau- Ginsburg theory using the same singular 

gauge transformation ||3^, Alternatively, one can also apply the statistical transmutation in order to keep the 
fermionic character by attaching an even number of flux quanta, which results in a Chern-Simons problem of interacting 
electrons at integer filling factors [|8[ |3|] . These attempts are guided by the desire to derive the quite successful trial 
ground state wavefunctions for the FQHE filling factors from first principles. Furthermore, field theories open the 
opportunity to calculate various observables in a systematic manner. The fermionic transformation was also applied 
to the case of even fractions v = l/2q for integer q by Halperin et al. ||2^. In this case, attachment of 2q flux quanta 
per electron, which point oppositely to the external magnetic field, leaves the system in average in zero magnetic field. 
This observation is the starting point for a reasoning that the system behaves effectively similar to a Fermi liquid in 
a vanishing field . An attempts to get rid of some inconsistencies in these theories led recently to the development 
of a Hamiltonian formulation theory by Shankar and Murthy [^l| . 

D. Theoretical models of the quantum Hall regime and the lowest Landau level approximation 

Our discussion hitherto emphasized that in the region of filling factors, where the FQHE occurs, interaction as well 
as disorder are necessary ingredients of any trustworthy theoretical attempt to understand the transport properties. 
On the other hand, the concomitant treatment of disorder and interaction is one of the great challenges that solid- 
state physics encounters. In fact, the difficulties are considerable as the example of the treatment of interaction in 
the problem of the metal-insulator transition in disordered solids shows Much worse, in most cases each single 
problem on its own is not exactly solvable. 

In the FQHE regime, the interaction is much larger than the fluctuation of the random potential since otherwise 
the many-particle gap above the ground state would be washed out and the effect would disappear. Because of the 
cleanness of the samples, disorder should not play the decisive role as long as transport theories are not considered. 
Nonetheless, disorder will quantitatively influence the results. In any case, our starting point is an interacting model 
without disorder. At strong fields, moreover, an independent particle model is an inappropriate starting point at partial 
filling of LLs due to the macroscopic ground state degeneracy. This argument does not question the usual treatment of 
transport in the IQHE regime, where one starts from a one-particle model with disorder neglecting just the interaction. 
The situation at integer filling factors is different from that at partially filled Landau levels as the unperturbed ground 
state is non-degenerate. Nevertheless, any consistent theory at integer v has also to account for the influence of the 
interaction, which has the same order of magnitude as the cyclotron energy under typical experimental conditions as 
we will see below. There are some recent attempts to incorporate electron-electron interaction in connection with the 
study of the diverging localization length in the center of the LLL indicating a localization-dclocalization transition 
in a random potential . Eventually, this whole discussion gets an additional interesting turn due to approaches 

that relate the system at a partially filled LLL v — p/ q to a model of weakly interacting new composite particles at 
higher integer filling factor p. 

Based on these assumptions, we can write down a simplified second quantized Hamiltonian describing interacting 
electrons, which move in a plane and are subject to a perpendicular magnetic field B = Be^, i? > 0, 



H = 




(8) 



The first line renders the one-particle Hamiltonian of Eq. (pj) using an effective mass m and gyromagnetic factor g of 
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the host semiconductor. The total spin density operator 



(T.CT' 



can be expressed by means of the field operators ^<T(r), ^'^(r) for electrons with spin a and the Pauh matrices f; [i is 
the chemical potential, and N = J^a I '^^''^i(r)^'o-(r) is the particle number operator. The second line describes the 
interaction between the electrons, where A = e^/(47re^c) is the coupling constant of a Coulomb like interaction. The 
interaction term takes into account the two-particle interaction V{r — r ') between the electrons in the plane as well 
as the interaction of these electrons with a homogeneous, compensating neutralizing background charge rif, and the 
self-interaction of the background. The sum of these terms can be expressed by means of the normal order symbol 
: and the electron density operator n(r) = ^J.(r)^o-(r). The charge neutrality is satisfied via the constraint 
J (Prn{r) = Uf,. In the following, we assume an isotropic particle interaction V{\r — r'|). 
Despite the restriction to the plane, the interaction is three-dimensional as the lowest subband wavefunction (f>o{z) 
is extended in the z-direction. Hence, the Fourier transformation of the interaction V^(q) is corrected by a form factor 
when the influence of the finite width w is discussed. The constant subband kinetic contribution to the Hamiltonian 
is omitted. 

The purely electronic Hamiltonian Eq. (^) still covers a huge number of physical situations ranging from the two- 
dimensional jellium model for vanishing magnetic field | |4^ to the strong magnetic field situation, where the QHEs 
occur. If we neglect the Zeeman energy for sake of simplicity, the ground state properties in the thermodynamic limit 
without magnetic field are characterized by the dimensionless electron gas parameter = rg/aB-, where rg = l/y^Tm 
is the average particle distance in two dimensions and gb = ineh^ / (me^) is the Bohr radius. If we switch on a 
magnetic field, the ratio of Coulomb to cyclotron energy 

'^^p^ = (10) 

determines the ground state properties. In the case of an infinitely large magnetic field leading to a vanishing 
interaction energy against the cyclotron energy, it is Vs — > 0. Then, the ground state properties become only a 
function of which is just the situation we encounter in the LLL approximation. Low values of mean high 
densities n. Consequently, strong magnetic fields are needed to reach the low filling factor region. Note that for a 
typical density of n ~ 10^^ cm~-^, it is ~ 1.8. 

Since ic/^o ~ \J v jl^ it is appropriate to consider the coupling constant A for filling factors ~ 1 as the relevant 
interaction energy scale. If one is interested in the thermodynamic properties of the model, the following three energy 
scales of an ideal 2DES govern the physics. These are the Landau level distance fiuj^ = h\eB\/m, the Zeeman energy 
Az = \gfJ,BB\ — \gBe\Ti/ {2m) and the Coulomb coupling constant A = e^/(47re£c)- For GaAs, the three energy scales 
expressed as temperatures in K are 



huJc/kB [K] = 20.05 B[T] 
A^/fcs [K] ^ 0.2953 B[T] 
X/kB [K] = 51.44 /B[r] , (11) 



when the magnetic field strength B is given in T. For magnetic fields greater than 6.58 T, the Landau level distance 
exceeds the interaction energy scale, whereas the latter one is still larger than the Zeeman energy up to an exper- 
imentally non-accessible huge magnetic field of 30 344 T. Although most of the strong magnetic field experiments 
are done in the range between 7 T and 20 T and at temperatures below 10 K, the omission of higher Landau level 
mixing is not automatically justified if the filling factor is smaller than two. Nowadays, e. g., samples with small 
carrier densities down to n ~ 0.2 x 10^"'^ cm~^ are available, and therefore QH-experiments &t v = 1 can be performed 
in a magnetic field of 2.5 T Then, the Coulomb energy exceeds the cyclotron gap, and the LLL approximation 
has non-negligible corrections. However, we restrict ourselves to the lowest orbital Landau level throughout this 
work, since the strong magnetic field effects become more pronounced at stronger fields. Furthermore, there are some 
technical simplifications as the form of the wavefunctions in the LLL that suggest such an approximation. 

The Hamiltonian Eq. (||) can be restricted for < < 2 to the Hilbert space of the LLL orbital states n — 0, when 
the condition 

A^,A<nwc (12) 
holds, and when the experimentally relevant temperature range satisfies the inequality 

fcsT < hujc ■ (13) 
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In first quantization, the Hamiltonian reads without a positive background term in the case of spinless electrons 
(0<:^< 1) @ 

H=^N +\j (0^>(q) « - n^=,e'^'''^) , (14) 

where riq — ^ ^'^'^^ ^^e Fourier transformation of the one-particle density operator, n„ = e"*'^^^^ e"'' 
means its projection onto the LLL, and z = x — iy,q — qx — iqy are complex coordinates p2| , 

Furthermore, we quote two second quantized representations in planar geometry. The first is suited for the treatment 
of a finite system on a disk in front of a positively charged background 

TT14=0; cr,cr' 

T7i,m'=o, m,m'— 

Here, N ^ ^ Ni is the particle operator, and the field operators are represented solely by the LLL eigenfunctions 
(pm{'^) and LLL fcrmionic field operators c^.o-, cln ^- ^■ 

00 

^'^W = "P^riir) Cra,a ■ (16) 

m=0 

The interaction matrix element in symmetric gauge for the LLL is given by 



Vr™„„„™3,„, = J d^J d'r'^*^^{r)^*^^{r')V{\r-r'\)^^,ir')^^Ar) (17) 
= {mi,m2\V\m4,m3) 

and satisfies for any isotropic interaction potential the angular momentum conservation mi +m2 = + 7714. Usually, 
we will omit the constant kinetic energy hujc/'^ multiplied by the particle number. In Eq. (p^), the constant neutralizing 
homogeneous background density nt is represented by the expectation value of the filled orbital LLL weighted by the 
filling factor v, which ensures charge neutrality and for models with long-range interaction the finiteness of quantities 
in the thermodynamic limit. 

In Landau gauge, the LLL Hamiltonian of Eq. (jlj) neglecting the kinetic energy term reads 

H = -Ia,{N^-N^) + ^ y Wiq,p-p')c\, c^ , ,c, , ,c , , (18) 

2 2 P -5:0- p +5,0- p~f,(T ^ ' 

where W(q,p — p') is the Landau gauge matrix element 

W{q,P-p') = I ^^(k)e-^e''=^(P-''')^^4„,, 

^ + (19) 

using the one-particle eigenfunctions of the LLL in Landau gauge 



o,fc,,(r)-^==-e (20) 



with centers Xi = 2-Kill/Ly {I = 0,±1,±2, . . .) along the x-direction. By excluding the q — contribution from the 
sum we take into account a positively charged homogeneous background. Diagrammatically, this corresponds to the 
neglect of the tadpole diagrams. 

Due to the macroscopic degeneracy of the non- interacting part of the Hamiltonian (p^, zero temperature field 
theoretical methods cannot be applied except for the special case i' — 1 Therefore, only a thermodynamic 
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Green's function method like the Matsubara technique allows a high-temperature expansion in X/{kBT) at arbitrary 
filling factors. For later use let us introduce the Matsubara or imaginary time Green's function defined as 

g,,,,{x,x')^-{T^^{x)^l{x')) , (21) 

where x = (r,T), T is the time ordering symbol, \E'^^^(a;) = e^^^^^'^/'' ^'^'''^(r) e"^^"^^'^/'' and (...) means averaging 
with respect to the grand canonical ensemble. The range of the imaginary time t is < r < h(3. For time independent 
Hamiltonians the GF becomes a function of r — r' and is antiperiodic with respect to the period h(3 in case of fermions. 
In a constant magnetic field the GF is independent of the degenerate quantum number k, i. e. 

^«,fc,<T;«',fc',<T'(^'^') = -('^C„^fc^<,(T)c,\, ,^, ,^,(t')) 

= 5n,n'6k,k'Sa,cr'Qn,a{T ~ T ) . (22) 

The free electron frequency dependent Matsubara function Qn}cr{'ivn') is given by 

1 1 

Si^ii^n') = = - / dre-"-eW (^) (23) 

with fermionic Matsubara frequencies Vn' — (2n' + l)7r/(?i/3) and electronic energies S,n}a = hujc{n + 1/2) ^ — /j, 
in the LL n with spin a. 

In case of the projection onto the LLL, the interaction energy scale dominates over any other energy scale except for 
the temperature, which makes the free electron problem inappropriate as a starting point for a perturbation theory. 
This problem is typical for the treatment of strongly correlated systems. Within the LLL approximation, the constant 
kinetic energy term with the band mass m can be neglected. Nevertheless, there exist situations, where the omission 
is not justified. For example, the construction of the current operator makes the inclusion of higher LLs necessary. 
Note that the LLL projection leads to an effectively one-dimensional problem. The Hamiltonian shows also some 
resemblance with lattice models of itinerant electron systems. The kinetic energy term can be viewed as the vanishing 
limit of a narrow band, but it will attain a non-zero value due to existing weak disorder. 

The LLL Hamiltonians (p^), (^8|) can be considered as a special case of a Hamiltonian describing a double-layer 
system consisting of two parallel quantum wells Such a system is characterized by the distance of the layers 

d ~ lOnm, which has the order of ic and hence of tq, and a hopping matrix element Ag^^s for interwell tunneling. 
For sufficiently narrow wells and strong fields, we can again assume that the electrons occupy only the lowest of the 
subbands in each well and equally only the LLL. Moreover, it is assumed that the electrons are spin-polarized. Then, 
in the unphysical limit of vanishing layer distance d = 0, but non- vanishing AgyiS: the system can be exactly mapped 
onto the one-layer system in the LLL with spin-degree of freedom and Zeeman term. The energy difference Ag^s of 
the symmetric and antisymmetric combination of the subband wavefunctions is the analog of the Zeeman energy A^ . 
This mapping is accomplished by a unitary transformation in the layer-index space introducing a fictitious pseudo- 
spin ||5l|, In case of a non-zero distance d > 0, the 5f/(2)-symmetry of the Coulomb interaction is broken to a 
t/(l)-symmetry as the inter-layer Coulomb interaction becomes smaller than the intra-layer interaction. This leads to 
a variety of zero and finite temperature phase transitions, which, for example, were experimentally and theoretically 
studied at total filling factor v = 1 |5^ . 

In this work, we are mostly concerned with ground state properties of the B ^ oo limit, where A <C A2, fuvc at 
the typical FQHE filling factor = 1/3 as well as with ground states and thermodynamics at and near v = 1. In 
the latter case, the experimentally important condition A^ <C A < TiUc holds. Therefore, we will ignore any effects 
related to higher LLs. We focus in the work on the bulk properties of the 2DES and do not touch two other important 
problems. First, due to the edges of the system there exist in contrast to the bulk region gapless excitations, whose 
low-energy excitations can be described as a chiral Luttinger liquid ||5^, |5^ . Second, one can discuss the physics of a 
system with a finite number of particles, which are confined to an approximately parabolic potential of small extent. 
This corresponds to the situation of a so-called quantum dot, whose physics is intensively studied without and in a 
magnetic field. 

Although this work considers only the physics at single points along the line of continuous values of the filling factor 
V, we will get an impression of the abundance of phenomena occurring in this small filling factor range Q < v < 2. 

Research in the field of 2DES in a strong magnetic field has been a lively and still expanding field for more than 
twenty years. Besides the original papers, there are review articles and books giving an introduction for the newcomer 
as well as valuable information for the active researcher. A still very useful source of information on 2DES up to the 
advent of the quantum Hall effects is the review by Ando, Fowler, and Stern from 1982 Q. Meanwhile, there are 
a few textbooks and collections on the integer and the fractional quantum Hall effect available ^ |63| . 
Pedagogical introductions to the field are lectures given by MacDonald Esf and Girvin iQ . 
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E. Two interacting particles in a magnetic field 



The study of two interacting spin-polarized fermions of equal and opposite charge, respectively, qi — ±(72 in a- 
magnetic field is very instructive for the understanding of the many-particle problem. 

The problem is further facilitated when the Hamiltonian is projected onto the LL n. Here we focus on the projection 
onto the LLL n — 0. But even pro jec tion onto orbital LLs n > is quite common when discussing the FQHE at, e. g., 
ly = 5/2 |65) or higher LL effects On the other hand, LL mixing cannot be ignored as was shown for 2 < < 4 
p7[. The two-particle Hamiltonian can be written as 




10.0 15.0 20.0 25.0 
Angular momentum M 



FIG. 3: The interaction energy spectrum of two spin-polarized fermions of equal charge in the LLL interacting via Coulomb 
interaction. The energy levels are drawn vs. the total angular momentum M — m_ + m+. Different energy levels indicate 
eigenfunctions with different relative angular momentum m_ whose eigenvalues are just the pseudopotential coefficients Vml 
characterizing the Coulomb interaction. Due to the antisymmetry of the entire wavefunction and the assumed spin-polarization 
the relative angular momentum eigenfunctions exhibit only odd m_ values. The horizontal degeneracy is due to the degeneracy 
with respect to the center of mass angular momentum eigenvalue. As long as the total angular momentum M < nimax, where 
rrimax was chosen to truncate the Hilbert space of the LLL, the breaking of the translational symmetry due to the finite 
extension of the sample cannot be read off from the spectrum. 

H = ffo.l + -ffo, 2 + -ffint, 1-2 

= Ho,i + Ho,2 ±XJ ■0^V{q)e^^^^^'^^'> , (24) 

where the sign depends on the repulsion or attraction of the particles' charges. V{q) is the Fourier transformation 
of the two-particle potential and equals, e. g., 27r/|q| for the bare Coulomb interaction in two dimensions. The 
separation into a free center of mass motion and the relative motion reduces the problem to a one-particle problem 
in the relative coordinate. In the case of equal charges and an isotropic interaction, the energetic eigenvalues after 
projection onto the Hilbert space of the nth Landau level are characterized by n and the relative angular momentum 
TO-. The eigenvalues are given by 

= hu;,{2n + l) + XVj;:^] 

= nu;,i2n + l) + X J ^^(g)e-|ql^^^ [L„(i^)]2L,„_ (Iqf.p) 

= hu;,{2n + 1) + «J„_ |i/™tKL J , (25) 

where the Ln{z) are Laguerre polynomials and the vJi^} are called pseudopotential coefRcients characterizing the 
projection of any isotropic interaction onto the Landau level n [ |68| . The eigenfunctions of the free relative coordinate 
problem in symmetric gauge are given by ip'^ln_ ■ For the LLL n = the explicit energy eigenvalues for the Coulomb 
interaction are 

= nu;, + AKn = hcj, + A (26) 
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using the notation Vrr, 



Vr 



(0) 



These pseudopotential coefBcients of the Coulomb interaction form a series of mono- 



tonically decreasing numbers, c/. also Sect. i. e. Vq — -s/tt/^, Vi = V^/^^, ^2 — 3-\/7r/16, V3 — 5^/tt/32, . . .. They 
can be visuahzed by numerically diagonalizing the LLL two-particle Hamiltonian exploiting the conservation of the 
total angular momentum M = to_ -|-to+, see Fig. ^. In contrast to the discrete energy values in the repulsive case, for 
attractive charges we get a continuous spectrum from the interaction term in Eq. ( p4[) when projecting onto the LL n. 
This is due to the fact that the exponential in can be separated in a part containing only harmonic oscillator like 
Landau level ladder operators and a part that depends on the total magnetic translation operator for two particles T. 
However, the latter is the projection of the generator of the translation operator P of the center of mass coordinate 
without magnetic field onto an arbitrary LL n. Since the translation operator has a continuous spectrum described 
by the two-dimensional vector k the energy spectrum in the LL n reads 

£;(")(k) = ?...(2n+l)-A J ^nq)e-^[L„(^)]^e-(-x>^)^? 

= ?iwc(2n+ 1) - Aa(")(k) , (27) 
where a'-"^(k) describes the binding energy of a magneto-exciton in the LL n. In the LLL, this energy becomes 



a(k) = a("nk) 



(2^)2 

27 



y(q)e- 



2 e 



(28) 



The second line shows that for ki^ > 1 the quantity a(k) describes the interaction between two well separated single 




FIG. 4: The binding energy a(k) Eq. (^9|) of two fermions of opposite charges in the LLL for the unscreened Coulomb interaction. 
Asymptotically, it is a{k) ~ 



particles of distance For an unscreened Coulomb interaction (no tilde) the isotropic binding energy is depicted 
in Fig. and equals 



«(k) 



-e * H—r) ■ 



(29) 



Here, Iq{z) is a modified Bessel function. The value q(0) = ■\/7r/2 can be identified as the exchange energy as 
comparison with the value of the Fock diagram shows [pO 




II. GROUND STATES AND QUASIPARTICLES OF SPIN-POLARIZED ELECTRONS IN THE 

LOWEST LANDAU LEVEL 



A. The Laughlin function as a quantum liquid ground state at filling factor 1/q (g — odd) 



This Section is mainly concerned with the ground state and its charged elementary excitations, the so-called 
quasiparticles (QPs) at filling factor v = 1/q. We restrict ourselves to the strongest FQHE filling factor v = 1/3. We 
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review the ground state properties of the LaughHn wavefunction and develop the notion of the quasiparticles. We 
illustrate these results on the basis of a generic model for the electron-electron interaction, the so-called hard-core 
model, which we will frequently use in modified form throughout this paper. In particular, we study the properties of 
the quasielectrons and discuss the merits of various explicit trial wavefunctions proposed for the quasielectron. Two 
theories that try to explain the occurrence of "weaker" FQHEs at filling factors v — p/q with p > 1 are discussed: 
the older hierarchical theory and Jain's composite fermion theory. We focus on the disk geometry because direct 
comparison with the Laughlin theory originally formulated in the symmetric gauge is possible and experiments are 
done in a planar geometry. On the other hand, results for finite systems on a disk are affected by the edge and show 
large finite-size corrections. These are disadvantages when one is interested in bulk properties in the thermodynamic 
limit, since numerical calculations in the spherical geometry don't suffer from edge effects. However, we think that it 
is worth to elucidate the problem in the disk geometry. 

Let us start with the main theoretical arguments supporting the correctness of the Laughlin theory for interacting 
electrons in a strong magnetic field at v = 1/ q. We assume in this Section that A ^ A^, TilOc , i- e. the electrons are 



restricted to the subspace of spin-polarized electrons in the LLL. This corresponds to the B ^ oo limit in Eq. (15) and 
allows the omission of the first two terms. Such an assumption is not a priori justified when comparing with typical 
experimental field strengths in En. ( [ll| ) since the inequality ( |T^ ) is not satisfied. But the discussion of the spin-degree 



of freedom following in Section III will show that such an assumption is justified, e. g., at v — 1/3. Therefore, in such 
a limit the only energy scale of the Hamiltonian remains the interaction strength A. 

The theory of the FQHE put forward by Laughlin p^, [7l]| is based on proposals for the ground state wavefunction 



of the interaction dependent part of the Hamiltonian (15|) at special filling factor v — 1/q and for the elementary 
charged excitations, the quasiparticles, at this v. At first glance, the search for an appropriate wavefunction seems 
to be hopeless as we encounter here a non-perturbative situation, where standard methods fail. However, Laughlin 
succeeded in finding a ground state wavefunction at the FQHE filling factor v ^ l/q [q - odd) with some remarkable 
properties. It can be shown that 

1. the Laughlin wavefunction describes a many-particle quantum liquid- like state, 

2. the elementary charged excitations have a fractional charge of magnitude ±e/g, and the sum of their energies 
defining the quasiparticle gap A^e-g/i remains positive in the thermodynamic limit. 

3. The resulting incompressibility ensures at exactly v ~ 1/q the position of the chemical potential in the many- 
particle gap above the ground state, which is a necessary condition for the peculiar transport properties. 

Laughlin's wavefunction for a finite number N of particles a,t v = 1/q approximately describes the ground state for a 
large class of model interactions. It is given by 

N ^ 

vi/i/,(zi, . . . , z^) = n - ^-^^ ''''' • (30) 



The density corresponds to a circular droplet about the origin of radius R ~ ^2N^ic- The Laughlin wavefunction 
has total angular momentum MN{q) ~ qN{N ~ l)/2. The unnormalized wavefunction is written in first quantization 
and ensures antisymmetry for fermions due to the occurrence of only odd q. The wavefunction is built out of the 
symmetric-gauge single-particle states of the LLL n — 0, c/. Eq. (^). Therefore, only powers of Zi, but not z*, occur in 
the non-exponential part of (|30|). Note that analogous wavefunctions can be found for the Landau gauge in the torus 
geometry [^| and for electrons moving in the field of a magnetic monopole on the surface of a sphere p8| . For a finite 
particle number, the filling factor definition v — N/Nq, for the thermodynamic limit Eq. (||) needs to be corrected by 



a constant K{i/) forming the "magic table" in jG^, so that 

= v-'^N + K{iy) (31) 

for any Laughlin filling factor i/ — 1/q. Starting with q — 1 and K{h' = 1) = 0, we get for a new filling factor 
(i/ ')-! = -I- 2 a correction K{iy') = K{v) + 2. Thus, it is K{v = 1/g) = -(?_- 1), so that N^=qN -{q- 1) is 
the corrected finite particle number definition. In fact, applying this scheme to (pO[), we have the following finite-size 
filling factor definition for the Laughlin wavefunction 

v={N- 1)/{N^ - 1) = (TV - l)/7rw. ^{N- l)/{q{N - 1)) = 1/q . (32) 

Here, rrimax denotes the maximum one-particle angular momentum in '^i/q measuring the spatial extension of the 
electrons or, mathematically, the highest power in any of the Zi as the polynomial part of the Laughlin function is 



homogeneous. Simple power counting in (30) shows indeed that the wavefunction has v = 1/q. The relation ([3l|) is 
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not only useful in order to identify the FQHE states for finite particle numbers at 1/q, but also to identify states as 
particle-hole symmetry FQHE states or states at filling factor v — p/q. 

We see that if we fix, e. 5., the coordinates of the particles 2, . . . ,N, the particle 1 is very efficient in keeping the 
other particles apart due to the qth power of the relative coordinates (zi — Zi). This essentially lowers the interaction 
energy. However, this qualitative argument depends on the filling factor and cannot be universal, since we expect at 
very small filling factors v or large q that a Wigner crystal and not the Laughlin wavcfunction is the ground state. 

We note that for q = I and any N the Laughlin wavcfunction is exact and independent of the kind of interac- 
tion, since all single-particle states of the lowest orbital Landau level are occupied. Therefore, the only possible 
Slater determinant agrees for the non-exponential part with the Vandermonde determinant in the coordinates Zi 
<I>i(zi, . . . , zjv) = ni>j"(-^i ~ ■^j) ^^'^ '^^^ written either in the form of Eq. ( pO| ) or in second quantization as 

1*1, TV) - \N-l,N-2,...,l,0)^c^^_,c^^_^...c\cl |0) . (33) 

Additionally, in the case of two particles requiring m^ax — q aX v ^ 1/g, one is inevitably led to ^' 1/^(21,^2) — 
{zi — ^2)^6"^'^^' )/4 independently of the interaction. Rewriting this wavcfunction as a product of a free particle 
wavcfunction in the center of mass coordinate with angular momentum — and a LLL wavcfunction e~'^~' 
with q — TO- shows the direct relation to the two-particle spectrum for n = in Eq. (p^). 

However, for q > I and arbitrary particle number this function remarkably withstands any analytical translation 
into the language of second quantization for larger N despite its apparent simplicity. This leads to difficulties to 
perform the thermodynamic limit N ^ 00, see e. g. and explains why to date a lot of our conclusions is based on 
finite-size numerical calculations. Indeed, the number of Slater determinants built out of one-particle states increases 
dramatically with increasing particle number TV. For TV = 10 and v = 1/3, e. g., there are 135 281 Slater determinants, 
whereas the dimension of the sub-Hilbert space of total angular momentum M = 135 is 246 448 . This emphasizes 
the strong correlation in the ground state at fractional v — 1/3 in contrast to the case — I, where only one Slater 
determinant, the Hartree-Fock state, occurs. 

Evidence that the Laughlin wavcfunction contains the essential physics comes from various sources. These are 
mainly 

1. the fact that the Laughlin wavcfunction is the exact solution of a model, the so-called hard-core model, with a 



specific interaction |68, |74| [75| 



2. the comparison of the numerically calculated exact ground state wavefunctions for a small number of particles 
"^ex with the Laughlin function ^'i/g in various geometries. Here, for small number of particles the energy 
expectation values were compared and the overlap 

ll*e.|| ||*1/,|| ^ ' 

was studied for Coulomb interaction and logarithmic interaction with an overlap larger than 0.99 fl^ , 

3. the analogy with a well studied classical problem, the so-called plasma analogy [ psf , 

4. field-theoretical derivations of the Laughlin wavcfunction after a singular gauge transformation, see below, by 
studying fluctuations either about a bosonic mean-field theory in an average zero magnetic field [^6|. |37|, or 
about a fermionic mean-field theory with an average magnetic field of the filled lowest Landau level ||3^, . Note 
that these field theories are justified by their capability to derive the already known ground state wavefunctions. 

The first argument pro vides some firm ground one can build on. Based on our considerations of the two-particle 



problem in Subsect. IE, the many-particle Hamiltonian for an isotropic two-particle interaction can be written in first 
quantization as 

N 00 

H,nt=XYl E ' (35) 

i>j = l 1=0, 

where the pseudopotential coefficients Vm"^ are given by Eq. (|25|), while the quantity \n,m)ij{n,m\ij acts as a projec- 
tion operator onto states with relative angular momenta to in Landau level n between particles i and j. Projection 
onto the LLL in ( |35| ) yields the real space representation 

N 00 

ffi°](ri, ...,r^)^Xj2 ^^T. ^™V2™j2(r, - r,) , (36) 

i>j—l m— 
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where 5'^{r) is the two-dimensional delta function. 

Besides the Coulomb interaction, a whole class of Coulomb like interactions can be constructed by truncating the 
series of pseudopotential coefHcients Vm at relative angular momentum k so that the coefficients Vo,Vi, . . . ,Vk are 
positive, but of arbitrary strength. In particular, they can agree with the values one gets for the Coulomb interaction, 
while Vm = if m > k. As already mentioned, even pseudopotential coefficients do not c ontr ibute to the energy 
expectation value, but they become important when the spin is taken into account, see Sect. III. In general, we call 
this model the hard-core model or more precisely in dependence on the truncation the Vk model. It is remarkable 
that the Laughlin wavefunction '^i/q is the exact zero energy eigenfunction of the Vq-2 model {q odd) for any finite 

N when neglecting the neutralizing background, i. e. the relation {'^i/q\Hy\,^\'^i/q) — holds. In the disk geometry, 
this eigenstate is unique and can be identified in the many-particle spectrum from finite-size diagonalizations, see 
Fig. 1^ ||78| . The reason is that only relative angular momentum components not smaller than q occur in the Laughlin 
function and therefore do not contribute to the energy as Vm ~ for m > q. To put it in other words, if the 
wavefunction vanishes at least like {zi — ZjY in each of the relative coordinates of two particles, the total energy is 
zero in the Vq-2 model. Reversely, if a relative angular momentum lower or equal to {q — 2) occurs in the wavefunction, 
which is the case at = l/g for all many-particle eigenfunctions except for the Laughlin function, the energy becomes 
larger than zero due to the positive contributions from Vi, . . . ,Vq-2- Therefore, the Laughlin wavefunction is for a 
given q among all possible zero-energy wavefunctions of the Vq-2 model that with the smallest maximum one-particle 
angular momentum rrimax- Hence, for v < 1/q, zero energy eigenstates, which are necessarily the ground states, are 
possible too. 

In particular, if we choose at ^ = 1/3 as generic model the Vi model with Vi — and Vm — for m > 3, the 

two-particle interaction is 



(37) 



A whole class of eigenfunctions with zero energy at filling factors v < 1/q can be constructed by multiplying the 



0) 

c 

UJ 




60 65 70 75 80 85 90 

Angular momentum M 



FIG. 5: The low-lying energy eigenvalues for the V\ model at = 1/3 and N — 7. The non-degenerate ground state of zero 
energy at MN=7{q = 3) = 63 can be identified with the Laughlin state Eq. (pOl). From W^. 



Laughlin state with a homogeneous symmetric polynomial S{zi, . . . , zjv) so that the new wavefunctions become 



v<l/q{zi, . . . , Zat) = S{ZI, ZN)'^l/qizi, ...,Zn) 



(38) 



Reversely, any wavefunction with filling factor larger than 1/q has a positive energy because in any case a relative 
momentum {q — 2) occurs. 

In reality, we cannot expect an interaction of such a kind to occur. Already in the case of the Coulomb interaction, 
the Laughlin wavefunction is not the exact ground state. This is due to the occurrence of non-zero pair am plit udes of 
relative angular momenta m > 3 (m - odd) in the numerically determined ground state wavefunction, c/. |80[ . Since 
for the Coulomb interaction all pseudopotential coefficients Vi, 1^, V5, . . . > occur, they contribute to the energy 
expectation value, which is different from the Vi model. Nevertheless, one can argue that the coefficients Vq with 
q > 3 are weak perturbations of the short-range model introducing additional smaller energy scales besides the energy 
scale set by ViA. This should not lead to significant differences between the exact eigenfunction found numerically 
and the Laughlin function. Thus far, this argument was never analytically pursued, but it is strongly corroborated 
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by extensive numerical calculations for finite-size systems. To check this, the overlap Eq. (34) was studied in various 



geometries by many authors. We mention for Coulomb interaction an overlap of 0.979 for = 4 in the disk geometry 



|_8, 81 , an overlap of 0.99406 for nine particles on the sphere |8^, and of 0.9881 for six particles in a planar geometry 
with periodic boundary conditions . 

The robustness of the Laughlin function with respect to changes in the interaction potential was shown for some 
typical interactions. However, there is also strong numerical evidence that the Laughlin wavefunction is not always a 
good choice for the ground state. Such a result is due to Haldane and Rezayi |Q who separated the LLL interaction 
Hamiltonian ( |3^ ) in the Vi-pa.rt and the residual part containing the influence of the pseudopotential coefficients Vm>i, 
i. e. Hintil) — + jHv,„yi- The numerical study of the low-lying spectrum varying the parameter 7 in Hintil), 
which comprises for 7 = 1 the Coulomb interaction and for 7 = the Vi model, revealed different behavior for 7 < 1.25 
and 7 > 1.25. In the former case, the ground state has a large overlap with the Laughlin function and is gapped, 
whereas in the latter case, the gap disappears and the Laughlin state is only found in the excitation spectrum. This 
observation emphasizes in spite of the artificial character of such interactions that ground states of repelling electrons 
in the LLL are not so universal like, for example, the BCS ground state in the theory of superconductivity, where 
already an attractive interaction between electrons leads to a new ground state. In summary, the occurrence of 
the Laughlin state as ground state depends on parameters like the filling factor and the pseudopotential coefficients 
characterizing the interaction. These conditions can be satisfied in experiments at certain fractional v. 

Now, we turn to some of the remarkable properties of the Laughlin ground state, which can be found upon studying 
the two-point distribution function 

(2)1 l\ * I \ * I l\ I l\ I N (*l/9l"^nil'^ni2''™2'^"''J*l/9^ 
n^'{z,z)= 2^ (^„^(z)(^„^(z (Z (z) — 1^-^ 

™' ll^l/qll 



N^oo J d^Zl...J d^ZN\'^l/q(Zl, . . . ,ZNjr 



(39) 



and the one-particle density 

^^'Hz) = jT--, I dz'n^-'\z,z')= >' <^(z)v.,„,(z) 



i (fz2...<PzN\^l/q{z,Z2,...,ZN)\^ 

= 1™ ^ -n^ r-pi — TT — 7 rP7 (40) 

N^cx, j (Pzi... \ CpZN\'^l/q{zi,...,ZN)\^ 

for the unnormalized Laughlin wavefunction '^i/q- The two-point distribution function is related to the pair distri- 
bution function by g{z^ z') = n^'^^{z, z')/n^. 

At v = I, the exact two-point correlation function in the thermodynamic limit 



n 



^'Hz,z')^^{l~e-l^-^'\-/') (41) 



exhibits translational invariance and isotropy and results in the homogeneous one-particle density n^^\z) = 1/(2it). 
For q > 1, the homogeneity of the Laughlin wavefunction, i. e. n'^^^(z) = i'/{2tt), can be inferred from a relation of 
^i/q to a well investigated model of classical statistical mechanics. By rewriting the squared modulus of the Laughlin 
function as the partition function we encounter the Hamilton function of the two-dimensional one-component plasma 
(2D0CP) jT8|. The plasma parameter F, which is proportional to the inverse temperature of the classical problem, 
can be identified with F = 2*7 As for q < 70 the 2D0C0P describes a fluid, while for g > 70 a crystal occurs, this 
yields a rough estimate for the transition between a fluid-like ground state with constant density and the crystalline 
behavior identified with the occurrence of the expected Wigner crystal. However, direct evaluation of ( ^0|) for q = 3 
up to 25 particles in the disk geometry shows quite large fluctuations around the expected constant density value 1/q 
pGf . In contrast, the ground state at magic filling factors in the spherical geometry has even for finite N a constant 
density distribution on the surface of the sphere due to its identification with a state of zero angular momentum L = 

The radial pair distribution function ^(r = |z — z'|) for odd q > 1 was analytically calculated by determining the 
coefficients of a power-series in |z — z'p using charge- neutrality, perfect-screening and the compressibility sum rules 
as constraints |87||. These results agree quite well with data from Monte Carlo (MC) calculations for q = 3, 5 up to 
144 particles |^8|, ^9{. The radial pair distribution function g{r) shows the correlation hole at r = and tends to one 
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for r ^ oo, which is typical for a fluid in contrast to the oscillating behavior of g(r) for a crystal exhibiting long-range 
spatial correlations. 

The knowledge of the distribution function and the radial pair-distribution function g{r), respectively, allows to 
calculate the ground state energy per particle e(iy) in the thermodynamic limit 

A r \u r°° 

^TTTr d?^' V{z, z') (z, z') - n') = — drr V{r) {g{r) - 1) . (42) 



2iV . 

At I' = 1, we find for the Coulomb interaction and the Vi model, respectively, 

econKl)/A = -^yi^ -0.6267, (43) 

evAl)/X = 2V, = ^. (44) 

Similar calculations at g > 1 show a lower energy of the Laughlin wavefunction when compared with other trial 
wavefunctions like, e. g. , charge density waves as long as g < 9 The calculation of the pair-distribution 

function of the 2D0CP at T = 6, i. e. q = 3, allows to determine eco«i(l/3)/A = -0.4100 ± 0.0001 which is 
almost identical with the value found from an MC evaluation of the Laughlin wavefunction in the disk geometry 
ecotii(l/3)/A = — 0.410 ± 0.001 Moreover, these results can be compared with finite particle number data from 
exact diagonalization studies at = 1/3 in the torus geometry with the extrapolation value —0.409510 from 

diagonalizations in the disk geometry j?^ , with the value —0.415 ± 0.005 Q in the spherical geometry, and 
with the result —0.4056 ||5^ based on the hypernetted-chain method of the 2D0CP. Similar results can be found for 
!/ = 1/5. Although these results show the superiority of the Laughlin wavefunction compared with other wavefunctions 
describing a crystalline state, this does not ultimately rule out other states with lower energy in the thermodynamic 
limit. 

Because of the restriction to the spin-polarized LLL, particlc-holc symmetry holds, which relates in the thermody- 
namic limit the ground state energies per particle at filling v and (I — v) by the relation 

(1 - 1^) e(l -~iy)- ve{v) = (1 - 2v) e(l) (45) 

regardless if a neutralizing background potential occurs. 
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FIG. 6: The ground state energy per particle vs. filling factor for the Coulomb model with neutralizing background in the disk 
geometry using data for six particles. Here, the filling factor definition v = N/N^ is employed. The calculated points are shown 
as filled circles and the long-dashed line serves as a guide to the eye. No precursor of a cusp is seen at = 1/3 for this finite 
particle number. However, extrapolations N ^ oo for 1/3 and two adjacent values might indicate the expected cusp when the 
extrapolated data (open circles) are connected by a solid line. Note the large finite-size corrections of ~ 13% at v — 1, where 
ecouiil) = -0.6267A. 
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B. Elementary excitations: quasiholes and quasielectrons as quasiparticles and the incompressibility of the 

Laughlin state 

The ground state properties of the Laughhn state are remarkable on its own. However, in order to explain the 
unique transport properties an additional ingredient is necessary, viz. the occurrence of gapped neutral particle-hole 
excitations above the ground state (l8[ |8|, |9^. This is equivalent to a discontinuity in de/dv at v — 1/q creating 
downward cusps in the e{i') dependence or to the incompressibility of the Laughlin state as we will see below. Cusps 
of e(v) at, e. g. V = 1/3, as they can be seen in finite size diagonalizations, c/. Figs. ^ and ^, are of limited value 
and do not prove the existence of a gap above the many-particle ground state. In order to show this, we employ 
the quasiparticle construction by taking properly the thermodynamic limit into account. The quasiparticles (QPs) at 
y = \/q are related to the ground state on the left and right side, respectively, infinitesimally close to 1/g, and they 
can be viewed as the elementary charged excitations of a quantum Hall system [Q. There are three ways to create 
such QPs by slightly changing the filling factor v = {N — l)/{BA/^o — l) = l/gtol/q±(5i/by changing either the 
area A (neutral QPs), the particle number N (gross QPs) or the magnetic field B (proper QPs). 

Within the construction of neutral QPs, the system area is either increased or decreased by an amount of 27r£^ 
altering the number of elementary flux quanta = BA/^q by one elementary flux quantum = h/e while keeping 
N and B fixed. The new ground state when compared with the homogeneous Laughlin state exhibits either a localized 
charge deficiency (one fiux quantum added) or a localized accumulation of charge (one flux quantum removed) of total 
charge \e\/q whose location is parameterized by the coordinate The former QP is called the neutral quasihole (QH, 



index — ) and the latter the neutral quasielectron (QE, index -|-).|20C] The total energy difference between the new 
ground state and the ground state at the QHE filling factor defines the neutral quasihole energy e"(i/) and the 
quasielectron energy e" (i^), respectively. These energies are given in a finite system and in the thermodynamic limit, 
respectively, by 

el(N,N^)=Eo(N,N^±l)-Eo(N,N^) 



e'»^T.M'|j^^. (46) 



There are two other ways to create QPs, viz. gross QPs with energy e±(v) and proper QPs with energy e±(v). In 
the first case, the number of particles is changed by one and leads to the formation of q independent gross QPs in 
order to conserve the total charge, while the latter case is based on the change of the magnetic field and therefore 
altering A'$ by ±1, but keeping N and A fixed. In each of the three constructions, the relative charge is increased or 
decreased One should note that these QPs are different from quasiparticles used in Landau's theory 

of Fermi liquids, where they are introduced to describe an interacting fermionic system approximately in terms of 
non-interacting new particles, whose life-time diverges like when approaching the Fermi surface ||95| . 

One can easily show that in the thermodynamic limit the various energies are related by the ground state energy 
per particle e(i/), g || 

el (v) = e± (v) T ^y) = l± (v) ± \ve(v) . (47) 

Furthermore, in case that the single QPs are identified in a QHE system, we can construct neutral excitations above 
the QHE ground state consisting of a QE and a QH far apart and hence non-interacting. The energy of such a 
particle-hole pair is given by the sum of their energies regardless of the kind of construction as Eq. ( |4^ ) shows 

A,e-..(.) =.;(.) + .'!(.) ^.^((I)^^-(I)^^ J . (48) 

Note that the positive definiteness of the gap implies downward cusps of the e(i^) curve as suggested in Fig. ^because 
of de/dv < 0. An instructive example of gross QPs can be found for non-interacting electrons at v = 1, see Fig. ^. 
The addition of one electron leads to a ground state with one electron promoted to the minority spin n — Landau 
level, i. e. the gross quasielectron energy is e+(l) = huJc/2 + A^. On the other hand, the removal of one electron 
from the ground state in Fig. |l| causes only a change of the kinetic energy so that the gross quasihole energy is 
e_(l) — —fujjc/'^- Therefore, the quasielectron-quasihole gap equals lS.qf.-qh(v = 1) = which is indeed the energy 
of the lowest-lying particle-hole excitation. Note that the non-interacting system is for non-integer v compressible 
since /S.qe-qh(v < 1) = 0. But this is altered at FQHE filling factors when interaction is taken into account. 

To be definite, let us first turn to the quasiholes of the Laughlin state illustrated by numerical results &i v — 1/3. 
In order to find an approximate analytical expression for the ground state of the QH wavefunction, Laughlin gave 
a prescription for the construction of such a state [Esl. The construction resembles the way a proper QH is created 
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as we increase the magnetic field in average by one flux quantum. Let us assume a one-particle energy eigenstate 
^m^-\z\ When a magnetic solenoid with flux $o and magnetic field parallel to the external magnetic field is 
added to the one-particle Hamiltonian, a new eigenstate 2;™+ie~l^l of same energy can be found after gauging 
away the vector potential outside the origin. If such an operation is applied to the non-degenerate Laughlin ground 
state in symmetric gauge "^i/q by adiabatically switching on a solenoid at the origin, the Laughlin state will evolve 

approximately in a new many-particle state ^[/^j where the angular momentum of each particle is increased by one. 
The resulting many-particle wavcfunction 



N 



,zn) 



J|zi^'i/^(zi, ...,zn) 



(49) 



is particularly efficient in keeping away the electrons from the origin and describes a quasihole at this location. 
Comparison of ( ^ ) with Eq. ( |38| ) shows that this quasihole is an exact eigenstate of zero energy of the Vq-2 model 
due to the symmetric prcfactor zi • . . . • zjv- Moreover, there exist in total N of such symmetric polynomials that 
increase within the neutral QH construction the maximum angular momentum by one, i. e. v — {N ~1)/N,^ < 1/q, and 
which are zero energy eigenstates of the Vq-2 model. This TV-fold degeneracy of QHs with total angular momentum 
reaching from M^r (g) -I- 1 = M* -I- 1 to Mjv {q) + N = M* + N and the Laughfin state itself, see Fig. |^, can be exploited 
to construct a QH located at ^ by superposing these -I- 1 states with weights ^™ (0 < to < N) 
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(zi, . . . ,ZAr) = Y[(^ 
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k>j = l 



(50) 



An alternate way to derive the QH wavcfunction at is the application of the magnetic translation operator (^) to the 
QH located at the origin (|49|). This leads to a shift of the QH at ^ = to the position ^ in front of the homogeneous 
background 



N 



N 



,Ziv) = *l7i^(zi,...,z^)e-^l«l^/^ne^ 



.r/2 



(51) 



Laughlin's QH wavcfunction can be recovered by neglecting the last exponential, which means projection to the 
correct filling factor y%. The charge deficiency is localized around ^ within the magnetic length scale what can 
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FIG. 7: The low- lying excitations in the V\ model showing the !^ = 1/3 QHs for a seven particle system. Besides the Laughlin 
state at Mff^riq = 3) = 63, each of the zero energy eigenstates is non-degenerate and is described by the states of Eq. ( |5C| ) 
with fixed total angular momentum M ranging from 64 to 70. From |79]. 



be analytically checked aX v — 1 since 



^W(z) = ^(l-e-l^-^lV^) 



(52) 



In total, it amounts for one charge, while far away from ^, the constant density of the Laughlin state remains 
unchanged. For v = 1/q, q > 1, the one-particle density is given by the expression ^2) multiplied by ly with a total 
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charge diminished by 1/g when comparing with the Laughlin state. The correctness of Laughhn's QH function for odd q 
and Coulomb interaction is again corroborated by generahzing the plasma analogy mentioned above. Now, additionally 
to the 20CP Hamiltonian function a term shows up, which indicates an interaction between the plasma charges q 
and a phantom charge of strength +1. Similar conclusions can be drawn from the overlap of the trial wavefunction 
with results from exact diagonalizations and from MC calculations. Exact diagonalizations are particularly successful 
in the spherical geometry because the QH eigenstates can be identified as degenerate states with angular momentum 
L = N/2 due to their {N + l)-fold degeneracy |6|. 

Another trial QH state formulated in second quantization takes the Laughlin construction literally [Q. Then, the 
operators [/,„ — c^„+iCm + 1 — rim are successively applied to the second quantized Laughlin state |^i/g) aX 1/q with 
as angular momentum distribution operator. It shifts each occupied one-particle angular momentum 
|m + 1). This corresponds to a QH at the origin and reads in second quantization 



state |m) 



) =C/o[/l...[/oo|*l/,) 



(53) 



It contains the same Slater determinants as the Laughlin QH (^), but this QH wavefunction differs in its weights 
of the Slater determinants. Despite the disadvantage that this QH state is not an exact eigenstate of the hard-core 
model, it is analytically more accessible and its energy expectation value for a more realistic interaction needs not 
to be energetically higher than that of Laughlin's QH state (pO|). This is indicated by Tab. ^ where for Coulomb 
interaction and Vi model QH energies and QE energies are listed at filling factor 1/3. 

When turning to the QEs, it seems at first glance that the concept of QH and QE creation is symmetric w. r. t. the 
stable Laughlin state at v = 1/q. However, there exist differences between these two types of QPs, which become 
already apparent when looking for the QE wavefunction in the Vq-2 model. Contrary to the QH, no exact solution for 
arbitrary particle number is known. Moreover, no wavefunction with rrimax = li^ — 1) — 1 can prevent the occurrence 
of relative angular momentum {q — 2), i. e. e^{l/q) > for any finite particle number. Hence, quite a lot of trial 
wavefunctions was proposed for the description of the QE emphasizing various aspects of the problem. 

Laughlin proposed a QE wavefunction centered at |]l8| , which is given by 
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(54) 



k>l=l 



The derivative operator is the adjoint of the product operator in (pO|). 

MacDonald and Girvin constructed a QE at the origin in analogy to the QH construction ( p3|) 

l^l}^ >=doo...rfl(l-^o)|*i/5 > . 



(55) 



1 ^ rim decreases all angular momenta by one except the one for m — 0, which 



Here, the operator dm = 
has to be projected out. 

A third t rial w avefunction is due to Jain | p6[ and is related to the composite fermion (CF) theory, which is discussed 
in Subsect. II C . Here, each of the total angular momentum components M* — N + m + 1 (m = — 1, 0, . . . , iV — 2) is 
given by 



'^l/q,M* -N+m 
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(56) 



where initially a wavefunction is constructed that contains one state from the LL n = 1 before the LLL projection 
operator Vq is apphed [ p7[ . 

A fourth QE wavefunction, which transfers the idea of the magnetic translation (^l]) to the QE at the origin, is the 
disk equivalent of a function originally proposed by Haldane for the spherical geometry It describes a QE 
centered at ^ 



JV 



(57) 
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For completeness, we list a proposal due to Morf and Halperin for a QE at the origin motivated by the plasma analogy 
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{Zl ~ Zj)(z2 - Zj) ^''^ 



(58) 



where A denotes the antisymmctrization operator. All wavefunctions describe a local accumulation of charge |e|/g 
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FIG. 8: Numerically calculated expectation values for Laughlin's (projected to the correct v\ our, and Jain's quasielectron 
trial wave functions arranged with increasing M in comparison with exact energies for A'^ = 9 particles in the hard-core model. 
Low lying states with total angular momentum M > M* are attributed to edge excitations. From |79l M. 



around the position ^, although differences in charge distribution in the core region of the QE exist ||8j, ^ . A formal 
pro of for the existence of fractional charges in the QH and QE wavefunctions can be performed via a Berry phase 
|P9| calculation based on the explicit form of the QP trial wavefunctions. Then, the Q P ch arge e* = ±|e|/g can be 
inferred from the phase the QP gains when adiabatically moving around a closed loop [100|. 

For small particle numbers, exact diagonalization data can be used to check the quality of the trial wavefunctions. 
Moreover, trustworth y tria l wav efunctions can be evaluated by MC calculations for larger TV not accessible to diago- 
nalization studies 
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|10J. Fig. g| shows for nine particles such a comparison assuming a V\ interaction in the 
disk geometry. While the approximate energetic degeneracy of the numerically calculated low-lying states allow their 
identification with QE states of angular momentum M in the range M* ~ N to M* — 1, the quality of Laughlin's 
QE wavefunction (|5J), Jain's QE ( p6| ) and our QE proposal (|57| ) become worse when moving from the origin to the 
edge of the system. The calculations show that among the proposed QE wavefunctions after extrapolation to the 
thermodynamic limit Jain's QE wavefunctions come closest to the finite-size diagonalization results found in the disk 
geometry | |9^ . A similar outcome was reported for the spherical geometry, where Jain's QE is energetically superior 
to Haldane's QE wavefunction [103 . Eventually, it is worth to note that for electrons on a disk with Coulomb inter- 
action the QE degeneracy could not be observed up to nine electrons probably due to strong finite-size effects, 
whereas on the sphere the QEs appear as ground states with angular momentum L = N/2 like the QHs and hence 
exhibiting an {N + l)-fold degeneracy The essential quantitative result is that the quasiparticle gap Age_g^ at 
ly — 1/3 has a value of approximately O.IA ~ 5y^B[T] K for Coulomb interaction. This is considerably smaller than 
the quasiparticle gap HlUc for non- interacting electrons. Similar calculations yielding smaller quasiparticle gaps were 
performed at i/ = 1/5, 1/7, which are in accordance with the observation that the FQHE at 1/3 is the strongest one. 
These theoretical results for Aq^-qh allow also comparison with the experimentally determined activation gap A.{v). It 
is found from measurements of the temperature dependent longitudinal resistivity p^x at s ufficiently low temperatures 
ksT ^ A in the plateau region of a FQHE filling factor using high-mobility samples llOSj . The activated behavior of 
the longitudinal resistivity is governed by the relation 



Pxx{T) 



rxx^ 



(59) 



where p^^^ is a constant, which is thought to be temperature independent. Assuming a spin-polarized ground state 
and non-spin-reversed QPs, which is the case at high enough magnetic fields, and which is supported by a non-zero 
Zeeman term, we identify A with A^e-g/i and hence a dependence Aqe-g/i(l/3) ~ O.IA oc is expected to occur. 
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Model, geometry 
method 


e 




4 ^- 






CouL, sphere, diag. |^ 


-0.410 


0.0947 


0.0089 0.2314 


-0.1278 


0.1036 


Coul., disk, analyt. ^ 
Eqs. (PI), (PI) 


-0.4098 


0.0959 


0.017 0.2337 
±0.002 


-0.120 
±0.002 


0.114 
±0.002 


CouL, disk, MC Hi 
Eqs. (poD, (pf 


-0.410 
±0.001 


0.0943 


0.0046 0.231 


-0.132 


0.099 
±0.009 


Coul., sphere, VMC [104] 

Eqs. (poD, (pep 










0.106(3) 


Vi model, disk, 
diag. PI 








0.1865 


0.1865 


0.1865 


Vi model, sphere, 
diag. PI 








0.1905 


0.1905 


0.1905 



TABLE I: The ground state energy per particle e, the neutral quasihole and quasielectron energies el and e^, the gross 
quasihole and quasielectron energies e_ and e+ and the quasielectron-quasihole gap Aq^-qh in units of X a,t v = 1/3. The 
results encompass data from exact diagonalizations as well as from energy expectation values of trial wavefunctions quoted 
by their numbers in the text (QH wavefunction first) for two models of spi n-po larized electrons in the LLL. Moreover, the 
geometry of the system used and the method applied are listed. In one case |104|, only the gap energy Aq^-qh is given, since 
the energy of an excited state with a QE and QH far apart was calculated. 



In experiment, this can roughly be observed, however, the gap value is at least more than 50 % smaller than the 
calculated ideal value at 10 T [105|. Effects reducing the activation gap are the finite extension of the wavefunction 
in the z-direction, disorder and Landau level mixing as we will discuss later in more detail. 

Before closing, let us me ntion two remarkable properties of the QPs. First, the QPs are anyons, i. e. they obey 
fractional particle statistics [106|. If in three dimensions two quantum-mechanical particles are interchanged, i. e. one 
particle arounds the other one in a half circle, the two-particle wavefunction acquires a prefactor of e*'^ — ~1 and 
gi27r — fQj- fermions and bosons, respectively. However, in two dimensions this phase can have continuous values ip 
with < (/? < 27r for topological reasons. To evaluate the phase for the case of two QHs, one employs the two-quasihole 
wavefunction Hil-^j ~ ~ £,2)^i/q and calculates the extra Berry phase due to the adiabatic movement of the 

QP located at |i around the other one fixed at ^2- This results for Laughlin's QP wavefunctions (|^, (jsj) in a phase 
factor for interchange oi ttv = tt / q 100 1 indicating that the QPs are n/q anyons. Note that the QHs of the v = 1 
state are fermions. Unfortunately, this property has no direct experimental consequences. 

The second surprising property is the already mentioned fractional charge ±e/q' of the \/ q QPs whose experimental 
detection caused some debate in conjunction with Laughlin's gauge argument |]l5| , 107, 108|. The most direct evidence 
for this fractional charge was found recently in two different experiments. In the first one, resonant tunneling through 
an anti-dot was performed, where the observed periodic conductance peaks in dependence on magnetic field and gate 
voltage allow to determine the charge of the quasiparticles of the 1/3 state ]109[ . In the second experiment, a non- 
equilibrium property was exploited, viz. the shot noise of the backscattered tunneling current through a constriction 
between two v — 1/3 quantum liquids. Th e shot no ise is proportional to the charge of the carriers and experiments 



showed indeed an effective charge of e/3 [110, 111|. In order to exclude the possibility that such an experiment 



measures rather the condu ctiv ity but not the QP charge, a similar experiment was performed at v — 2/5. Here, a QP 
charge of |e|/5 was found |112|, which is in accord with the prediction that QPs at v = p/q carry charge ie/g. 

Now let us turn to the incompressibility of the Laughlin state. The isothermal compressibility can be either micro- 
scopically calculated as the long- wavelength limit of the density-density correlation function or from its thermodynamic 
definition 



1 



(60) 



At zero temperature, this can be related to the ground state energy per particle t{v) via 
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Incompressibility is therefore achieved by a discontinuity of the first derivative of the ground state energy at a certain 
filling factor, cf. (61), which corresponds to the expected cusp in an e{v) representation or by a jump of the chemical 
potential /i, which differs in its left and right values iJ,± = {dE/dN)_^_ = {d(ye(y)) / dv) j^. In a finite particle system, 
this jump in the chemical potential A/i = /i+ — /i_ is 

l:\li{N)=E{N + 1)+E{N -l)~2E{N) , (62) 

however, the thermodynamic limit is decisive as our discussion of the cusp in Fig. ^ in a finite particle system showed. 
From the comparison with (48), it follows that the jump of the chemical potential is proportional to the quasiparticle 
gap 
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(63) 



Therefore, the existence of a positive quasiparticle gap Aq^-qh means incompressibility. 

In contrast to the compressible gas of non-interacting electrons without magnetic field, the situation changes for 
non-interacting electrons in a constant magnetic field when a Landau level is filled. For arbitrary j/, the ground state 
energy per particle in units of hujc depends only on ly, even for a finite particle number p9| because of 



e(^) 



1 



(64) 



At integer filling u — n, the jump of the chemical potential, which is located in the charge gap in the middle of the 
two adjacent LLs, becomes A/j, = 1/n, and therefore it is Aqe_g/j = 1. Thus the system is incompressible for each 
integer while it is compressible for any non-integer v. Note that the quasiparticle gap is independent of the Landau 
level and particle number and can also be found from the finite N expression (62). Of course, the reason for the 



incompressibility is the existence of discrete Landau levels, and this holds also when the electron's spin is taken into 
account. Now let us return to the case of interacting electrons, where we expect also for FQHE v incompressibility. 
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FIG. 9: The ground state energy per particle evi vs. filling factor v for the Vi model without background. For < < 1/3, 
the energy equals zero (dashed curve). At 2/3 < < 1, we find evi(i') = (2 — l/!^)evi(l) (dashed curve) from ji^), where 
evi(l) = 2ViA is the ground state energy for the filled LLL. In the range 1/3 < < 1, the crosses are the result of finite-size 
diagonahzations for seven particles on a disk indicating strong finite-size corrections when compared with the exact results in 
the range 2/3 < i' < 1. In order to iniprove the results, extrapolations have to be performed for each u. The self-consistent 



Hartree-Fock approximation, see Sect. [II, yields a linear dependence evi iy) = 2yV\\ (0 < < 1) and shows its limited worth 
except at u — 1. Note that the quasielectron energy cannot be estimated from the slope at = 1/3 in contrast to equivalent 
data from the spherical geometry. The energy is given in units of e^i (1) = 2ViA. 

We showed in the last Subsection that this is the case at ~ 1/q by extrapolating diagonalization data and trial 
wavefunction results to the thermodynamic limit. In the Vi model at v — 1/3, such a quasiparticle gap develops only 
by the non-vanishing contribution from the quasielectron energy. This is illustrated in Fig. ^. On the other hand, this 
intuitive argument for a finite particle system cannot replace a formal proof since the probability amplitude of this 
contribution remains unknown in the thermodynamic limit. Up to now, a rigorous analytical proof hasn't been found 



yet 1 13, When comparing this situation with that for non-interacting electrons, we realize that we know neither 
the exact QE wavefunction nor a lower bound of its energy. There are other attempts to show the incompressibility 
at I' = 1/3 in the Vi model. They are based on a diagrammatic high-temperature expansion of the grand potential 



113|, a direct high-temperature expansion of the energy |114 , and a virial expansion [115 
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C. The microscopic understanding of filling factors p/q 



Besides the strong FQHE aX v = 1/g for g = 3,5 and their particle-hole analogs at z/ = 1 — 1/g in the LLL 
approximation there exist much more v, where plateaus in p^^y emerge, cf. Fig. ^j. The most prominent ones are of the 
form V = p/{2p+l) < 1/2 and the corresponding particle- hole symmetric states v = l—p/{2p+l) = {p+l)/{2p+l) > 
1/2 for p — 2,3,4. So far, we outlined the main features of the Laughlin ground state strongly corroborating the 
existence of an incompressible ground state at v = 1/q. However, this does not explain the existence of ground states 
with similar properties at v = p/q for p > 1 like the 2/5-state. 



The hierarchical theory |25[ |68|, 80, |89|, ^ tries to apply the stability condition at = 1/q (|32) between the 
number of degenerate electrons N and the degree of degeneracy to the equally degenerate boson-like quasiparticles 
of the parent state at 1/q. If the quasiparticle filling fraction attains the value 1/p {p - even integer), the 
hierarchical theory predicts two new second level hierarchy FQHE filling factors for the daughter states, either at 
V = pI {qp+ 1) = l/{q + l/p) in case of quasiholes or at = p/ {qp— 1) = l/(q — 1/p) in the case of quasielectrons. In a 
short notation these finite chain fractions correspond to a [g, +p] and a [q, —p] state, respectively. The denominators 
remain odd since p is even. Within this theory, e. g., the 2/7 state is the daughter state of the quasiholes of the 1/3 
state, and the 2/5 state is created by the condensation of the 1/3 quasielectrons if p = 2. 

In principle, this mechanism can be applied again and again by constructing new states from the quasiparticles of 
one hierarchy to get a stable state of the next hierarchy level unless the new incompressible state becomes unstable 
against the formation of a different kind of state as the Wigner crystal in the low filling factor region. In any case, 
the states in the hierarchical theory have filling factors with odd denominators. 



Clearly, there are some shortcomings in this theory |116|. First, the starting point of the construction are states 



at 1/q, what seems to be rather arbitrary. Second, the number of quasiparticles becomes macroscopically large in 
order to reach a new filling factor. The deeper one goes in the hierarchy by consecutive application of the hierarchical 
construction the more questionable the results become. Third, the essential prediction of the hierarchical theory is 
that a daughter state is not observed unless all of its parent states are observed. This agrees with experiment and 
shows that the quasiparticle gap Age_g/j of the daughter state is smaller than that of the parent states. On the other 
hand, the experimentally observable state at v = 3/7 appears in the short notation as the [3,-2,-2] state, while 
I' ~ 5/13 corresponds to [3, —2, -1-2] and hasn't been yet observed. Naively, one would expect that both states have the 
same degree of stability. Fourth, there are fundamental difficulties to write down explicit wavefunctions. There were 
attempts to introduce wavefunctions using quasiparticle coordinates Josf , but there remains the difhculty to translate 
them into electronic wavefunctions. This problem is closely related to the occurrence of anyons and their fractional 
statistics, which the quasiparticles obey Last, but not least, there is also a technical difficulty omitted so far: the 
QPs are not point-like and therefore corrections are needed in order to render correctly their short-range interaction 



117| . 

Up to the general acceptance of an alternative and more accessible theory by Jain |115], there were, nevertheless 



few attempts to put the hierarchical theory on a more quantitative footing by determining the interaction of the 



quasiparticles on a sphere from finite size numerical calculations for the case of quasielectrons [117 and for the 



quasiholes |11£|. In the first case, the authors succeeded in predicting quite accurate values for the ground state 
energy and the quasiparticle gap by determining the QE interaction of the 1/3 parent state, whose QEs condense into 
the 2/5 state. In the latter case, the authors found that the low- lying excitations of a bosonic system with Coulomb 
interaction at quasihole filling factor v = 1/2 are similar to those of an electronic system at v = 2/7, whereas the 
ground state energy was not determined. In another attempt we considered the quasiholes of the v ~ 1 state in the 
Vi model and showed that the ground state energy of the electronic 2/3 st ate is identical to the ground state of an 
effective bosonic Hamiltonian of interacting quasiholes at filling factor 1/2 [|l2C| . 

In summary, on one hand, the hierarchical theory is at first glance qualitatively appealing, however, it is rather 
clumsy and of limited value in their predictability at general filling factors, in particular, when compared with the 
composite fermion theory. 

The composite fermion theory by Jain [|9^, |ll§| , |l2l| proposes explicit electronic trial wavefunctions for the ground 
state of N interacting particles at filling factor v = p/q {q - odd). The construction of the wavefunctions at a filling 
factor of the form = p/{2kp -I- 1), fc - integer, is done in the following way: accommodate all N particles in the 
lowest p (spin-polarized) Landau levels. This results in the wavefunction <I>p, which is the Slater determinant of 
non-interacting particles at v = p. Then, $p is multiplied by the (2fc)th power of the Vandermonde determinant $1 
of N particles. Contributions from higher Landau levels, that inevitably occur for p > 1, are projected out by the 
application of the projection operator Vq leaving us with the wavefunction 

N 

Xp/(2fcp+i)(^i, ■■■,zn)= Va<l>fl>p = ■Pon(^» ^ zjf'^^p . (65) 
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Actually, the number of flux quanta per particle or the inverse filling factor is l/v = 1/p + 2k = {1 + 2kp)/p. This 
is the so-called D-operation within the formal classification of operations to generate wavefunctions at certain FQHE 
filling factors starting from integer v = p. In order to generate wavefunctions at other i/, Jain introduced additionally 
the C-operation allowing to create a wavefunction aX 1 — v when the wavefunction at v is known because of the 
particle-hole symmetry. Eventually, the L-operation creates states at 1 -I- from FQHE states at v. For example, the 
2/7 st ate appears as the result of the D-operation from 2/3, which again is the result of the C-operation from the 1/3 
state p2^ . 

These trial wavefunctions are meanwhile well accepted since extensive numerical calculations show that their overlap 
with the numerically determined ground state wavefunctions is very large. Moreover, the ground state energy per 
particle found from trial wavefunction expectation values and diagonalization results show an astonishing agreement, 
e. g. &t V = 2/5 for ten particles, the Jain function yields —0.4693, whereas the exact value is —0.4694 [123|. 



Furthermore, QP wavefunctions can be constructed by creating first the QPs of filled LLs and then performing the 
same operations as described above. Such a QE wavefunction was already discussed in the disk geometry, where 



projection is even necessary for the QEs of the Laughlin state, cf. Eq. (|56|). The important outcome of all these 
calculations at various v = pj q with p > 1 is similar to that of the Laughlin wavefunction: the neutral excitation of 
a QH-QE pair is gapped, and it follows the incompressibility at these filling factors. 

It is interesting to note that even without projection onto the LL L on ly a small portion of spectral weight resides 



in higher LLs so that the excess due to kinetic energy is quite small [124 . Two properties are particularly appealing: 
first, the Laughlin wavefunctions appear as a special case by starting from the filled LLL, i. e. , p = 1, fc = 1. Second, 
the unique wavefunctions for integer and fractional v are related by a simple operation, and a deeper connection 
between the FQHE and IQHE becomes obvious below. 

As in the case of Laughlin's function, the ultimate reason for the incompressibility is poorly understood. While for 
V = \ jq the hard-core model with the Laughlin state as exact solution gives a hint how incompressibility is caused, 
aX V = pI (2kp -|- 1) a similar model, where the Jain function appears as the exact ground state, is missing. However, 
there exists for filling v = p/{2p + 1) an interesting mechanism that supports the incompressibility property of the 
Jain function. In case of p — 2, we choose the somewhat artificial model, where the Hilbert space consists of only the 
two lowest Landau levels n — 0,1 and whose energetic distance vanishes, i. e. this corresponds to the band mass limit 
m ^ CXI. It turns out that in this case the function X2/5 is an exact solution of the Vi model, which is also related 
to a wavefunction for electrons with spin [ p5[ . It was shown by numerical calculation that in the limit m — > 0, which 
eventually leads to the LLL approximation, the overlap of the Jain function and the numerical diagonalization does 



not vanish and becomes approximately 0.73 |125| 



One rather puzzling feature, that we will discuss below, is the fact that the starting wavefunction was found for a 
Hamiltonian, whose energy scale is set by hujc, while the resulting Jain function describes the physics of a Hamiltonian 
of the interaction energy scale A. 

Now let us approximately reverse the above construction in a more physical way. We saw that the Jain function 
Xp/{2pk+i) provides an extremely accurate description of the ground state correlations of the interacting Hamiltonian. 



The statistical transmutation mentioned in Subsect. [C is in fact a singular gauge transformation with a parameter 
a counting the number of fiux quanta, which are attached to each electron and corresponds to an infinitely thin 
magnetic solenoid. A usual non-singular gauge transformation leads to an additional phase factor e*^^/'''-^ = e^^^^/'^° 
in the one-particle wavefunction when the vector potential A is changed to AA = A -I- V/ = A -I- a. In a singular 
gauge transformation, an additional magnetic field b = V x a occurs only at the location of the particle. In a 
many-particle formulation, such a transformation can be accomplished by adding, e. g. for the particle with index 1, 
a vector potential 

where for arbitrary i the condition i < j prevents double counting. The Qij = arctg{{yi — yj)/{xi — Xj)) is the 
azimuthal angle between Zi and Zj relative to the x-axis. The solenoid occurs for Zi = Zj with fiux a^o, where the 
sign of a determines the direction of the additional magnetic field relative to B. Note that (^) is a constraint between 
the vector potential and the particle density that is responsible for the Chern-Simons term occurring in a Lagrangian 
formulation of the problem ^ . Applying such a singular gauge transformation to the Jain wavefunction with the 
special choice a = ~2k yields a new electronic wavefunction obeying the same Hamiltonian except for the additional 
vector potentials a(zi). It reads 

Xp{zi, . . . , zpf) = Yl e "'^''^"Xp/i2kp+i) = n ( _ /.I ) Xp/(2fcp+i) • (67) 

i>j = l i>j ^ ' * ^ ' ^ 
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This leaves us in average with a magnetic field B* of strength B* = B — 2kn^o, where B = n^o/v = n^o{2kp+l)/p = 
B*{2kp+ 1), and n is the particle density. The resulting integer filling factor v* — p oi an inhomogeneous field satisfies 
the relation 

(68) 



2ki^* + 1 



We ignore the possibility that the resulting field B* is antiparallel to B, which would give ;/ = v* /{2kv* — 1) in (p3q). 

It is obvious that the gauge transformation is not exactly the inverse operation as can already be seen in the 
Laughlin case p = 1: we get neither the wavefunction of the filled LLL at i/ = 1, since we leave this subspace, nor is the 
magnetic field homogeneous. But this transformation is a starting point for field theoretical considerations to derive 
the Laughlin wavefunction [^s] , ^ by carrying out a gauge transformation at = 1/g with an even number of (g — 1) 
fiux quanta. The resulting fermionic field theory is evaluated at = 1 within a random-phase approximation (RPA). 
Also a bosonic formulation is possible, when a = —q is odd [^ . This is equivalent to the statistical transmutation 
applied to the Laughlin wavefunction that was used to prove algebraic off-diagonal long-range order (ODLRO) [Q. 
In the bosonic field theory, the particles move in an average zero magnetic field, and an approximate derivation of 
the Laughlin state is also possible. Moreoever, the field theoretical formulation opens the opportunity to calculate 
correlation functions. 

Now let us return to Jain's wavefunction, whose product Ili<j ~ ^j)^''*^p interpreted as weakly interacting 

"quasiparticles" in Landau level p, to which 2k fiux quanta are attached despite the difference to the singular gauge 
transformation (|67[). Indeed, the prefactors behave like vortices, whose arounding yields just a phase 2k x 27r. Hence, 
the idea arouse to consider the complicated many particle state of electrons at magnetic field B as weakly interacting 
composite particles, the so-called composite fermions (CF), in a reduced field B*. In this picture, most of the 
interaction is compensated by the 2fcth power of the Vandermonde prefactor $i. However, one should keep in mind 
that the term composite fermion is not only used in connection with Jain's construction, but also when Chern-Simons 
theories are actually meant. 

The idea of the Chern-Simons theory can be extended to the so-called "unquantized Hall effect" filling factors 



like 1/2. The fermionic field theory by Halperin, Lee and Read [g2j and their refinements, cf. |40, 126|, start by 
applying the gauge transformation with a — —2, which keeps the fermionic character, and consider the interacting 
particles in a vector potential fluctuating about zero. The main difficulty is to take appropriately into account the 
electron-electron interaction as well as the fluctuations of the gauge fleld, which couples to the electron density pO[ . 
The appealing feature of such a theory is the fact that the unperturbed ground state describes independent fermions 
forming a Fermi sea, while the original problem has no perturbational starting point. This theory predicts a Fermi 
surface and compressibility despite the perturbations. Hence, the transport behavior is very different from that at 
QHE filling factors. When the Chern-Simons idea is properly translated to Jain's construction of wavefunctions, we 
are led to a = 1/2 trial wave function of the form 

Xi/2^Vo\FS)Y[(z,-z,f , (69) 
where \FS) describes the Fermi sea of spinless electrons at zero field. The overlap of this Rezayi-Read wavefun ction 



with diagonalization results is again impressive in the spherical geometry, e. g., for nine particles, it is 0.9988 [127|. 
Of course, the existence of a Fermi surface should have clear experimental consequences. Experiments at the primary 
series of filling factors ly = p/{2p ± 1), whose limit for p — > cxj is just ly = 1/2, are elucidating. According to Jain, 
the effective magnetic field B* is related to the field aX v — 1/2 by ±i3* = Bp/(2p±i) ~ -Bi/2 — ^n^o/p- Hence, the 
activation gap energy Aqe-qh should scale like 1/p or linearly with B* as long as the effective mass m* is independent 



of B* |128|. Moreover, the symmetric appearance of oscillations in p^x around v = 1/2, cf. Fig. ||, resembles those 



known from Shubnikov-de Haas oscillations and replicate the features due to the IQHE a,t f = p. This suggests to 
interpret these oscillations as those of composite fermions in the weaker field B* . Thus, the FQHE appears as the 



IQHE of composite fermions M, 118 



Other peculiarities of the conductivity axx at and near v = 1/2 can only be seen when going beyond d. c. transport 
measurements by exploring length scales that are much sma ller than the mean free path £* of the impurity scattered 



composite fermions, i. e. for q£* ^ 1. Therefore, Willett ||20|, 129 1 exploited the surface acoustic wave (SAW) technique 
to get information on axxiq,^) 

via the amplitude damping and the velocity shift Avs/v^ of the SAW at sufficiently high frequencies (v^ = uj/q - 
velocity of sound). These quantities can be related to the conductivity, which, on the other hand, shows for ql* ^ 1 a 
linear g-dependence. At sufficiently high frequency (or wavevector q), the semi-classical theory in the small effective 
field B* predicts minima in the velocity shift Avs/vg of the SAW, which occur at certain values of the product qR* 
p2|. Here, i?* is the composite fermion cyclotron radius i?* = vp/uJc = fikp/{\eB*\). Hence, the observation of the 
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principal resonance (minima) in the velocity shift away from the minimum at -B1/2 allows the determination of the 
composite fermion Fermi wavevector kp, which is found to be very close to the value kp — \/47m of a gas of free 
electrons with one spin direction This result dramatically supports the CF picture. The origin of the SAW 

resonances is of geometrical nature due to the commensurate resonance between the cyclotron radius R* and the 
SAW wavelength 27r/g. Another experimental set-up studies the longitudinal resistance Rxx near = 1/2 in a regular 
arrangement of anti-dots, which acts like a superlattice with lattice constant a. Geometrical resonances are observed 



if a commensurate condition between a and the effective cyclot ron r adius R* is satisfied |130|. These resonances are 



similar to those found in real low magnetic field experiments 131 . Despite all these successes, the Chern-Simons 
theory has to answer the fundamental question of the correct energy scale, which is influenced by B* as well as the 
CF mass m*. The CF energy gap h\eB*\/m* oc \B\/{{2kp+ l)m*) must agree with the FQHE activation energy 
gap for every p in the primary series, which scales in the strong field limit like ^/B. Thus, the effective mass m* has 
a ^/B dependency at fixed v and has to diverge when approaching i?i/2 as p — > 00. This is difficult to achieve in 
the diagrammatic treatment by Halperin et al. . Therefore, the effective mass m* was used as a phenomenological 
parameter, which is fitted to the known activation gaps at a certain v of the primary series. Furthermore, experiments 



trying to measure to* are ambiguous and depend on the method used to extract the effective mass |132| , 133 



On the other hand, the unrenormalized band mass to cannot be completely abandoned. For example, it couples 
to spatial inhomogeneities of the magnetic field and enters the magnetic moment. Therefore, any satisfying theory 
should distinguish these kinds of masses. 

An approach that tries to circumvent these difficulties was recently developed |^l], 134|. After performing the 



fermionic statistical transmutation, Murthy and Shankar enlarge the Hilbert space by imposing additional constraints 
in order to properly consider the physical degrees of freedom. Approximate treatment of this problem and back 
projection onto the physical space yield Jain's and the Rezayi-Read wavefunction for ;/ — p/(2p + 1) and v = 
1/2, respectively. To go beyond the derivation of wavefunctions, canonical transformations and the approximate 
treatment in the long-wavelength limit result in a mass renormalization to higher values to* and exhibit the desired 
\/B dependency of the LLL approximation. These features are very pleasing and pave the way to a more unified view 
of the odd and even denominator states in the QH-regime. 

III. GROUND STATES AND THERMODYNAMICS OF ELECTRONS WITH SPIN DEGREE OF 

FREEDOM IN THE LOWEST LANDAU LEVEL 

A. The spin degree of freedom 

In the last Section, our considerations were restricted to spin-polarized electrons. However, the underlying assump- 
tion of an infinitely strong magnetic field is not satisfied in experiment. In contrast, the Zeeman energy is much 
smaller than the interaction scale, but the magnetic field is so strong that the inequality ( p^ ) holds. Hence, we return 
in this Section to the Hamiltonian of the form ( [l5| ) and (|8|), respectively, taking explicitly into account the Zeeman 
term. Then, the appropriate parameter to describe the various physical regimes is the effective gyromagnetic factor 
g = AJX cx Vb. 

This fact, however, does not devaluate our previous investigations in the limit g 00. The results are correct as 
long as ground state and charged excitations remain spin-polarized, even when (|l2|) is not fulfilled. This is possible due 
to the influence of interaction on the spin-orientation. Indeed, finite-size diagonalizations indicate that the Laughlin 
ground state at v — 1/g is even in the absence of a Zeeman term spin-polarized in the LLL approximation, while 
the ground state at v = 2/5 yields a spin singlet at vanishing Zeeman term, but becomes spin-polarized under the 
influence of a Zeeman term at accessible values for g in GaAs . Therefore, the existence of spin polarization in the 
ground state must be checked in each case. Even more intriguing is the question whether both charged excitations or 
only one of them is accompanied by spin reversals. Ideally, this problem can be solved by studying the magnetic field 
dependence of the activation gap A in the longitudinal resistivity pxx ~ e^^'^^'^, cf. Eq. (p9|). Then, a square root 
dependence on B indicates purely interaction driven excitations, whereas a linear B-dependence shows that spin-flip 
excitations are involved. Such different regimes could be observed, i. e., at lower fields a linear dependence on B with 
a cross-over to a square root -B-dependence when increasing the strength of the field. Unfortunately, the infiuence of 
disorder and the finite extension of the 2DES in the z-direction spoil often the experimental results. 

Trial wavefunctions at fractional ly considering the spin degree of freedom were early proposed ||25|| . But only with 
the extension of the LLL Hilbert space up to 1/ = 2, the physics at = 1 appeared in a new light. While in the spin- 
polarized case, the QE at v = 1 does not exist, now, the QE has to reverse at least one spin. Surprisingly, the charged 
excitations are not necessarily one-spin flip quasiparticles changing only locally the homogeneous spin magnetization 
of the ferromagnetic ground state at v = 1. At a sufficiently small effective gyromagnetic factor g, charged spin 



texture excitations with a large number of flipped spins are energetically preferred |135]. These excitations are called 
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skyrmions. Their existence has experimental consequences for the spin magnetization so that a lot of experimental 
and theoretical activity was sparked. 

In this Section, we are mainly concerned with the physics at exactly and near ly ^ I. The methods we apply are in 
some sense specific to this filling factor. This is in accordance with the fact that there exists no many-particle theory 
in the quantum Hall regime comprising the physics of the entire filling factor range. Therefore, general attempts are 
not able to concomitantly reveal the specific features of the physics for all v. Theories trying to explain experiments 
over a large range of filling factors as, e. g., tunneling experiments between two pa ralle l lay ers contain so strong 



approximations that they do not incorporate features specific to FQHE filling factors [136, 137 1. 

Now, the particle-hole symmetry for spin-polarized electrons is generalized to particle-hole symmetry around ly ^ I, 
cf. Eq. (p5|). For the ground state energy in the thermodynamic limit, e{h'), it holds 

(2 - ly) e(2 - ly) - ly e{iy) = 2(1 - z/)e(2) , (70) 

which is invariant with respect to the substitution ly ^ {2 — v). Since we are mostly interested in magnetic properties 
due to the spin of the electron, let us quote some of the repeatedly used elementary relations. Assuming an external 
magnetic field in the z-direction, the temperature dependent extensive total spin magnetization M is proportional 
to the difference of the expectation values of the particle numbers for spin-up electrons and for spin-down 
electrons 

Af (T, v) = \^{N,- N,) = Afo^^^ - Mo(^T - ^i) , (71) 

where Mq = N\g\iiB/'^ is the total magnetization of N spin-polarized electrons, and Va- is the temperature dependent 
filling factor for the respective spin direction cr with ly^ + lyi = v. The expectation value of the local spin operator 
( |l6| ) can be expressed for a state \^) in the LLL by 

SW = ^ E ¥'™(r)¥'™'(r)(*|4,.T;_,'C™',.'|^) (72) 

(m denotes here a general degenerate quantum number in the LLL) and defines the dimensionless local magnetization 
m(r) ~ 2S(r). Consequently, the total spin magnetization is 

M = / SrmAv) = j d\S.{r) . (73) 

When discussing spin configurations, it is interesting to know the total number of reversed spins K in the minority 
spin band. This number can be extracted from mz(r) via the relation 

J d\il-m,{r))-^{N.^-N) . (74) 

The interacting Hamiltonian ( |l8|) exhibits a continuous S'?7(2)-symmetry in the spin-space when the Zeeman term 
vanishes. Hence, the role of any non-zero Zeeman coupling is that of a symmetry breaking field. In case that the spin 
of the ground state for a vanishing gyromagnetic factor g is different from S ~ 0, the system develops spontaneous 
magnetization, which is parallel to the direction of the external magnetic field. We will see that such a situation 
can occur for quite a lot of filling factors in the LLL. If one compares this phenomenon with ferromagnetism in a 
metal, it is at first glance astonishing that in a QH-system such a situation occurs although the electrons are exposed 
to a strong magnetic field when the Hilbert space is restricted to the LLL. However, it is important to realize that 
the role of the symmetry breaking field is played by the entire Zeeman term, which can be independently tuned 
by the magnetic field as well as by g. Thus, the symmetry breaking Zeeman term can even vanish in the strong 
magnetic field limit. This can be realized in experime nts u nder hydrostatic pressure. This allows to sweep the g 
factor from its "natural" negative value through zero [ |138[ . Under such circumstances, the strong magnetic field 
enters only indirectly through a "band" of vanishing bandwidth and the LLL interaction matrix elements, which 
are quite different from the unprojected Coulomb matrix elements. Hence, a QH-system that exhibits spontaneous 
magnetization is called a quantum Hall ferromagnet (QHFM). One should keep in mind that the two-dimensionality of 
the system and the continuous SU (2)-symmetry for a v anis hing Zeeman term prohibits any magnetization at non-zero 
temperature because of the Mermin- Wagner theorem ]139( | . Any reasonable theory of the spin magnetization should 
satisfy this constraint. 
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B. Ground state and charged excitations at and near = 1 



We start with the discussion of the ground state of interacting electrons with spin in the LLL for fiUing factors 
exactly at v = I. 

The physics at = 1 is comparably well understood. We choose the Vq model for the interaction and construct 
a many-particle electronic wavefunctio n of zero energy in the spirit of the derivation of the Laughlin wavefunction 
at v = 1/3 for the Vi pseudopotential |140|. Note that due to the additional spin degree of freedom denoted by x 
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FIG. 10: The spectrum aX v = 1 for six particles with spin interacting via the Vb model on a disk (Vo = \/ti/2). The degeneracy 
of each energy level can be read off from the spectrum due to the chosen non-zero value of the Zeeman energy = O.OIA. 
E. g., the zero energy ground state at AI — 15 is split into seven non-degenerate levels, i. e. S = 3. 

the entire wavefunction has still to be antisymmetric, but the symmetry condition for the orbital part is relaxed and 
allows now also relative angular momenta m that are even. We know that any iV-particle eigenfunction of the Vq 
model without neutralizing background and Zeeman term of the form 

N 

M^,x]^[lliz^-z,)]^B[z,x] , (75) 

i<j 

has zero energy because of the occurrence of relative angular momentum one in the antisymmetric prefactor. The 
residual wavefunction for N particles with spin 1/2 has to be symmetric. The maximum power of Zi in 

the prefactor is iV — 1 so that for i/ = 1, i. e. N = 7V$, no further orbital dependence of as a power of z may 
occur (except for the exponential part). In order to be an eigenfunction of S^, vE'b[x] must be a symmetric A^-particle 
spinor with total spin S = N/2. The zero energy level is {N + l)-fold degenerate with spin component Sz between 
—N/2, —N/2 -I- 1, ... , N/2. For Vb > 0, any state different from ( [75|) has a positive energy, and therefore 

|*.=i) = |TiT2T3...Tjv)|*s) (76) 

is a zero eigenenergy state with maximum Sz = N/2, \'^s) is the Slater determinant of the filled LLL. All the other 
degenerate ground states with spin-eigenvalue 15,52) can be constructed by subsequent application of the ladder 
operator = E^i -5^ to @) as [S',S^] = 0. Any non-zero Zeeman term lifts the degeneracy and selects from 
the spin multiplet the state Eq. ( [76| ) as the unique ground-state of energy —NAz/2, see Fig. |l^. In the ground-state, 
all available orbitals with to = 0, 1, ... , iV$ — 1 are occupied, and the spin for each orbital state is up t- Hence, the 
schematic representation of the ground state in Fig. |^ remains exact for interacting electrons in the Vb model. It 
is clear that any double occupancy of an orbital costs energy due to the occurrence of relative angular momentum 
zero. Actually, the maximum Sz ground state is identical with the Laughlin state for q = I in the spin-polarized 
model. There, however, such a state occurs due to the infinitely large Zeeman term, while in the current case this is a 
pure interaction effect resulting from the spontaneous breaking of the symmetry. The Vq model without neutralizing 
background at = 1 plays a role analogous to that of the Vi model for the spin-polarized case at ly ~ 1/3. Its real 
space representation is in accordance with ( ^6|) 



H,nt = 47rVbA^ (5^(r, - r^) , 

i<j 



(77) 
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where we adopt for Vq its Coulomb value in the LLL, i. e. Vq = \/7r/2- The Fourier transformation of { fn\ ) is 
V{q) = AttVo so that a(0} = 2Vo, cf. Eq. @. 

A translation of Eq. ( |7q ) into an equivalent formulation within a bosonic language for ^'^[zjx] is possible. Then, 

the bosonic operators bnla with angular momentum m and spin a allow to write the degenerate ground state as 



|noi =K,noT=N- K) 



1 



^{N - K)\K\ 



K 



N-K 



|0) 



(78) 



where K = Q, . . . is the number of reversed spins as Sz — N/2 — K, Uma is the occupation number for bosons 
with angular momentum m and spin cr, and |0) is the bosonic vacuum. The state with K = corresponds to (|7^). 
Hence, the bosonic ground state is completely characterized by the occupation numbers uq^ of orbital states with 
m — 0. The situation does not change qualitatively when we consider the more realistic Coulomb interaction with 
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FIG. 11: The spectrum aX v — 1 for ten particles with spin interacting via the Coulomb interaction on a disk with a neutrahzing 
background. The Zeeman energy is Az = 0.005A. The ground state has, as expected, total angular momentum M = 45 and 
5 = 5. 



pseudopotential coefficients monotonically decreasing from the maximum value Vq. Although we cannot rigorously 
infer that the ground state is identical to that of the Vq model, finite-size numerical calculations in different geometries 

. In physical terms, the existence of aligned spins is favored by the Pauli principle. Parallel 



show that, see Fig. 11 



spins exclude the occupation of the same orbital state and thus reduce the overlap between two electrons and the 
Coulomb energy contribution. Moreover, parallel spins lead to a gain in exchange energy and to a decrease of the 
ground state energy compared with that of non-interacting particles. The mechanism is particularly effective for a flat 
band, whereas in a system of finite band-width a competition between increase of kinetic energy and exchange energy 
gain occurs when all spin are aligned. However, even the zero band-width is not sufficient to prevent the breakdown of 
the complete spin - pola rization leading to spin-singlet ground states when there is no Zeeman coupling, e. g. &l v — 2/?, 



and 1/ = 2/5 [141, 142 . In the Coulomb model, the ground state energy is the same as for the spin-polarized case since 
only the orbital part contributes, and we recover eco«;(l)/A = — a(0)/2 = —y^ir/S, cf. Eq. (|4^). In order to introduce 
the diagrammatic many-particle approach, let us repeat the ground state energy determination by calculating the 
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FIG. 12: Proper self-energy diagram in the self-consistent Hartree-Fock approximation. The propagator line in this diagram 
must be determined self-consistently. This approximation leads to a frequency independent self-energy and hence to a Green's 
function whose spectral weight consists of a single 5— function. The first order tadpole diagram is absent because of the 
neutralizing background charge. 
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thermodynamic Green's function (GF) in a self-consistent Hartree-Fock approximation (SHF), cf. Eq. (23), 



The combination of the Dyson equation with the temperature and filhng factor dependent, but frequency independent 
self-energy approximated by the Fock diagram yields 

On the other hand, the general expression for the temperature dependent energy |46| for electrons in the LLL within 
the SHF is 



6(°) + is.(zi.„) 



Y.r^F{^r)ie'"^+^r)- (81) 



At 1/ — 1, where particle- hole symmetry — —^^^ holds, the chemical potential is ^ = — Aa(0)/2. The zero 
temperature SHF energies for the spin-directions are e^^ = —A2/2 — Aa(0) and e^^ — A2/2, respectively. Assuming 

= 0, we recover the exact result for E{T = 0, i/ = 1)/A^ = eco«;(l) = -A&(0)/2 = -^/VJSX of Eq. (g). As 
usual, the consideration of the HF-diagrams gives for a uniform system at zero temperature a correction to the non- 
interacting result, which is half of the exchange energy contribution in the LLL 6°° = — Aa(0). Therefore, the SHF 
yields the exact ground state energy. 

After the discussion of exactly filling factor one, we turn to the charged excitations, the quasiparticles. Similar to 
the case of the spin-polarized ground states at v = 1/q, we can ask what happens if we remove one electron from, 
or add one electron to the spin-polarized ground state at = 1. However, in contrast to the spin-polarized case, 
where adding an electron a.t v = 1 leads inevitably to the occupation of the second orbital Landau level n — 1, 
and removing one electron creates iV -I- 1 quasiholes, which become degenerate in the thermodynamic limit, now, the 
emerging picture is much richer since we remain in the LLL utilizing the spin degree of freedom. 

We use again the Vb model for a microscopic understanding of the quasiparticles. First, we discuss the quasihole 
charged excitations, which are also much more amenable to diagonalization in the disk geometry as the data are not 
spoiled by finite-size effects. Nevertheless, particle-hole symmetry about v = \ makes these considerations not specific 
to the quasiholes. The neutral quasiholes are created by finding the ground state for N particles moving in an area 
covered by iV$ — N + 1 flux quanta. It is obvious that the various quasihole energies in the Vb model equal zero, 
i. e. e" (1) = 0, and hence the ground state energy in the Vb model (A^ = 0) is 




0<v <l 
Kv <2 



because of Eq. (70), see Fig. p_ 



On the other hand, the neutral quasielectron energy becomes in the thermodynamic limit 

e'l(l) = lim Eo{N,N^ = N -1)-Eq{N,N^^N) , 



,!»(l^l -'(|i| ■ (82) 



which leads to e" (1) = 2VbA , see Fig. |T^. In an alternate way, the finite-size quasiparticle energies can be also 
calculated directly from the energy of the quasielectron with a spin reversed particle in the state with angular 
momentum m, (m = 0, . . . , iV$ — 1), i. e. I^'^"'^) = cj^|cj| . . . cjy_2||0). In first quantization, the quasielectron state 
at the origin (m — 0) reads [|9| 

'i'[+l^,{z,,...,ZN)=A\ n (^.-^,) UiT2-..Tiv) I e-^S"il-'l' (83) 



with A as antisymmetrization operator. The representation allows, on the one hand, to identify the quasielectron 
states from the numerical calculations, see Fig. despite their strong finite-size effects, on the other hand, we 
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FIG. 13: The ground state energy per particle evo (v) as a function of v in units of evg (2) = VqA for the Vb model covering the 
whole u range < v < 2. Note the right side derivative aX v = of value 2, which determines the QE energy in the model. 
In contrast to this example, for Coulomb interaction only single eiiergy values are exactly known. In the case of spin-polarized 
electrons, approximate interpolation formulae for < < 1 exist (l43]. 



find in the thermodynamic limit {N^ oo) that the quasielectrons are degenerate and have energy 2VoA. Similar 
calculations for the conventional quasiparticle states in the Coulomb model with background lead to results sum- 
marized in Table ^ Below, we shall see that these states are not the true elementary excitations for the Coulomb 
interaction, i. e. energetically lower excitations exist. Furthermore, it is interesting to note that the existence of a 
non-vanishing gap for particle-hole excitations far apart is necessary for the incompressibility of a state. However, 
this does not exclude the existence of gapless neutral excitations occurring as long-wavelength spin-waves above the 
polarized ground state at = 1 in the case of vanishing Zeeman energy, see Subsect. |III C . This is different from the 
situation at fractional ly for spin-polarized electrons in the LLL, where even in the long-wavelength limit fc ^ of 



the single-mode approximation a gap exists |32 . It is remarkable that the case ly = I is the only example to date. 



Model 6(1) 6'I(1) eljl) 6-(l) e+{l) A,,,,,, A,fc_„,fc 

Vo model 2Vu 2Vo 2Vo 2Vo 
Coulomb -lyf lyf lyf ^ 1 



TABLE II: Exact results for the ground state energy per particle e(l), for the neutral and gross quasihole energies e" (1) and 
e_(l), the neutral and gross quasielectron energies e+(l) and e+(l), the value of the quasielectron-quasihole gap Aqe-qh, and 
the value of the skyrmion-antiskyrmion gap Ask-ask, see below, in units of A for the Vo model and the Coulomb model with 
background at = 1. 

where the incompressibility caused by interaction and not by a one-particle effect can rigorously be shown in the 
thermodynamic limit. From a microscopic point of view, the phenomenological distinction between integer QHE with 
a one-particle gap and the fractional QHE exhibiting a many-particle gap loses its meaning at = 1 as in the LLL the 
one particle gap is small in comparison with the enhanced spin-exchange gap due to many-particle effects. This 
becomes particularly obvious in experiment, where the exertion of hydrostatic pressure causing a vanishing (^-factor 



does not lead to the destruction of the QHE [138 . The classification of the charged e xcitation s of the quasihole type 



at 1/ = 1 for the Vb model is explicitly possible and was done by MacDonald et al. 140 , 144]. The QHs appear as 



zero energy eigenstates in Fig. [15L The multiple degeneracy of most of the states for a given total angular momentum 



M and Sz in the right panel of Fig. 15 shows that not all states have S = N/2, i. e. not all states can be described 
as spin split polarized quasiholes known from the spin-polarized case discussion in Section ||. A gain, all zero energy 
functions comply with Eq. (|75|). Hence, only boson occupation numbers npo- and riio- can occur ensuring that the 
maximum power of Zi is fixed due to A'^ = — 1. So-called seed states can be uniquely identified in Fig. ^5|, which 
are eigenstates of the angular momentum, of Sz and even of in the thermodynamic limit. Hence, the application 
of the ladder operator on the state with no| = K, i. e. |no| = K, rig-f = 0, rii^ = 0) vanishes so that this state can 
be uniquely identified as an energy eigenstate with S = Sz = M — N/2 — K. Then all states of the spin-split diagram 
can be identified since the degeneracy with respect to the center of mass motion leads to a degeneracy of states with 
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25 30 

Angular momentum M 



FIG. 14: The quasielectron at filling factor 1^ = 1 for N = 8 and N<s, = 7 in a model with Vb-interaction and no background. 
The energy is given in units of VqX. The lowest lying state at M = 21 has 5 = 3 with 7 particles spin up and one spin down. 
It has the lowest finite-size correction when compared with the thermodynamic limit of the energy of a quasielectron in the 
Vo model, which is £+(1) — 2VoX- The energy correction is the smallest one if the quasielectron has its innermost position 
minimizing the infiuence of the edge of the system. 
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FIG. 15: All skyrmionic states for six particles and N<s> = 7 at Az = O.OIA in the Vo model. The right panel indicates the 
degeneracy of each of the skyrmio nic st ates even after ordering w. r. t. the angular momentum M and in the presence of a 
non- vanishing Zeeman term. From |144|. 



given Sz along the M axis. Due to the mapping between bosons and fermions, the seed states can be translated into 
fermionic states with K reversed spins [145 
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(84) 



where li e {0, . . . , N}\{mi^ . . . , ttik}, h+i > h and C{K) is the normalization constant. 

This somewhat idealized picture in the Vq model is distorted when the interaction is of the Coulomb type so that 
the energy becomes lowest if the number of reversed spins K increases. This behavior can be quantitatively described 
by the energy of the generalized quasiparticles with K reversed spin €{K)± = E{K,N = iV$ ± 1) — Eq{N = A^$), 
where E{K, N = ± 1) is the lowest energy state with K flipped spins. Therefore, e{K) is monotonically decreasing. 
In a sufficiently large finite system, this means S = 0, see Fig. While in a microscopic quasiparticle state the 
additional spin has to change its direction on the microscopic length scale £c, just the opposite happens for large K, 
where the spins turn smoothly in order to minimize the exchange energy. Thus, it is tempting to derive a classical 
field theory for the spin magnetization m(r) with the local constraint m^(r) = 1 whose order parameter field changes 
slowly on the magnetic length scale ic- The order parameter field m(r) describes the deviation of the magnetization 
from the ferromagnetic gro und s tate m(r ) — (0,0,1). A gradient expansion of the energy functional £'[m(r)] yields 
in leading order two terms |135, 145, 146 1. The ffist term corresponds to a classical ferromagnetic Heisenberg model 
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and the second to a topological term, see below, 

[m] = ^ j <fr (Vm)2 - 2p, j d?r {m ■ [d^m x dyXn]) , (85) 
where ps is the spin-stiffness. The next order of the gradient expansion is of the Hartree type 

E.H = ^ /dV / d\' . (86) 



It describes the direct interaction between the local charge densities p{t) and is given at v = \ hy 

p{y) = -!-ea&m(r) • [dam{Y) x dbm{v)] . (87) 

OTT 

The mathematical structure of (^) ensures that the integration of the density p(r) over the entire two-dimensional 
space yields the topological charge 

Q= ( d\p{v) , 



which is always an integer number. Therefore, p{v) is called a topological charge density [147|. Magnetic field 
configurations of different Q can be easily distinguished by mapping the magnetization field from the plane onto the 
surface of the Riemann sphere. The integer Q counts how often the local field wraps the surface of the sphere showing 
that configurations with different topolo gical charges cannot be reached by continuous deformations. 

It has been known for a long time |l48[ that local extrema of (^5|) besides the global minimum energy field 
configuration m(r) = (0, 0, 1) characterized by the topological invariant Q exist. Such a classical field configuration, 
called either a skyrmion (Q < 0) or an antiskyrmion {Q > 0), is shown in Fig. |l6|. The size of a skyrmion Isk is defined 
as the distance £sfc between the origin of the skyrmion, where the spin points down, (mz(0) = — 1) to a point, where 
the z-component of the local magnetization vanishes, i. e. mz{r = Isk) — 0. Here, we employ the radial symmetry so 
that the magnetization vector is in the x-y-plane. Note the symmetry of the skyrmions with respect to translations 
of the origin as well as to rotations about an axis along the z-direction through the origin. In case of a non-vanishing 
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FIG. 16: Illustration of the local magnetization m(r) of a scale invariant skyrmion with a local magnetization at the origin 
pointing oppositely to the direction of the external magnetic field. The spin gradually turns up when departing from the origin. 
The extension of the skyrmion is defined by the radius Xah, for which raWW vanishes. From [ ^ . 

Zeeman term, we have to add the expression 

EWm]^^ I d\[l-mWr)] . (89) 



The explicit solutions minimizing the functional (85) are known for arbitrary Q, and a parameterization for the |(5| = 1 
skyrmion located at the origin is 

m(r) ^ (4^, -p^, , (90) 
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which leads to the topological charge density 



1 



p{r) = - 



4A2 



^(r2+4A2^)2 



(91) 



and Q — —1 because of {p8\j. Note that the topological charge is independent of the scale parameter A^fc. 

This independence is equally true for the energy of an antiskyrmion with topological charge Q > or a skyrmion 
with Q < 0, and the energy is given by 



(92) 



If we take into account the Hartree-term (p6|), a skyrmion of maximum size is preferred, but the topological charge 
and energy of the skyrmions remain unchanged. Therefore, without a Zeeman term, a neutral excitation consisting of 
a skyrmion-antiskyrmion pair with Q = ±1 has an energy Agk-ask — St^Ps = Aa(0)/2 = y^Tr/SA in case of Coulomb 
interaction. Most remarkably, this is only half of the quasielectron-quasihole pair energy Aqe-qh = Aa(0), see Tab. ||. 
On the other hand, in the Vq model, there is not such an energy gain for the skyrmion-antiskyrmion pair since 
Ps = Vo/ (47r)A and therefore A^k-ask = Aqe-qh = 2VbA. This is just what we found from the discussion of the highly 
degenerate ground states on the disk. In this respect, the Vq model is only instructive to link the microscopic picture 
with the description in terms of a classical order parameter field. 

Since the topological charge density and the electrical charge density are identical at = 1, the topological charge 
has direct experimental consequences. This is different from a Heisenberg ferromagnet, where only spin, but no charge 
degrees of freedom exist. The skyrmions can be observed in transport experiments. For example, characteristic features 
of skyrmions in tunneling are predicted in the current-voltage-characteristic, when electronic tunneling between two 
parallel 2DES near filling factor = 1 is performed [149|. 

In general, the Zeeman term (|8^) cannot be neglected, and this causes competition between a smooth change of the 
local magnetization with K ^ oo for vanishing Zeeman term and a single spin-flip when the Zeeman term dominates. 
The ground state is a function of the effective gyromagnetic factor g = A^/A. The function e(if (5)) determines the 
number of flipped spins in the ground state as well as the gap for charged excitations. For increasing ^, the classical 
model becomes less appropriate and the correct excitations can only be found from variation al w avefunctions, from 
a Hartree-Fock variational Ansatz, and from exact diagonalizations with different accuracy [ 145 |. The appropriate 
generalization of the classi cal field-theory is a formulation that takes into account quantum fluctuations around the 
classical solution [150, 151 1. 

Moreover, there exists a critical effective gyromagnetic factor Qcrit^ above which the skyrmions have a higher en ergy 
than the quasiparticles. This value was found to be gcrit = 0.054 from diagonalization studies on the sphere |13f;] 
and variational wavefunction calculations |1145| . From a diagonalization study in the dis k ge ometry, we got the value 



0.041, which is presumably smaller due to the strong finite-size effects in this geometry [144 . In GaAs, the magnetic 
field dependent ratio is g = A^/A = 0.00574-\/B[rJ. In typical experiments, it is g < 0.054 as in Barret t's NMR - 
experiment, where g ~ 0.016 [152[. Theory predicts for this value of g a number of ~ 3 — 4 flipped spins [145, 153 |. 



Note that although K should be quantized, such a behavior has not been seen in experiment so far 

Besides these findings, the charge spin texture excitations have the following other important properties: 

1. The charge, as we have already seen, and statistics of skyrmions agree with those of the corresponding quasipar- 
ticles of the Laughlin wavefunction. The 1^ = 1 skyrmions are spin 1/2 fermions and carry charge ±e [150, 154 



2. 



The occurrence of skyrmions as the lowest charged excitations at larger, odd filling factors depends on the kind 
of interaction. For example, at v — 3, g — 0, and zero width uu, skyrmions are not energetically favorable 
with a Coulomb interaction [155[, but for a width larger than a critical value and a g below a certain critical 
gyromagnetic factor, skyrmions occur as charged excitations |156{. 



3. There exists another critical gyromagnetic factor g^ 



10 ^. Below this value, two skyrmions can form a 
bound state, but the smallness of g^rit disorder effects presumably prevent any experimental observation 



4. Skyrmions exist possibly also at the fractional filling factor v = 1/3, however, g is rather small [158 



These theoretical results are corroborated by measurements utilizing different experimental techniques. Ground state 
spin magnetization measurements about ^ = 1 were among the first experiments, which gave evidence of the existence 
of skyrmions [152[, see Fig. |l^. In a simple phenomenological model, the slope of the magnetization curve at v = 1 
can be related to the number of flipped spins K of a skyrmion {v — 1~) and an antiskyrmion {v = 1"*"), respectively, 



cf. Eq. (106). Any K > 1 leads to a drop of the spin magnetization that is stronger than that for simple quasiparticle 
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excitations, see Subsect. [II D. From the slope of the curves in Fig. |l^, an estimate of the number of flipped spins gives 
K ~ 3.6 ± 0.3 in good accord with theory. However, some uncertainty remains, in particular, the spin-polarization 
value at r = 1.55 K is probably smaller than the assumed complete polarization for T = at = 1 |15£|. Second, 




FIG. 17: Knight-shift measurement data from an NMR-experiment by Barrett et al. in a GaAs-AlxGa\-xAs heterostructure 
at T = 1.55 K. The Knight-shift is proportional to the spin-polarization of the 2DES, and the maximum is due to saturated 
ferromagnetism at exactly v = 1. T he p recipitous drop of the spin magnetization is attributed to antiskyrmions on the left and 
skyrmions on the right side. From |152]. 



-^^ for sufficiently small temperatures in the QHE plateau region. The concomitant 
= B±^ fixed allows a continuous sweep oi g = A^/A to larger values while staying at 



magneto-absorption experiments, see Subsect. [II F|, allow the determination of the spin-polarization as a function of 
v. Values of K between 2.5 and 3.7 were foundpl. Third, as explained earlier, the standard method to determine 
the activation gap A for polarized FQHE states is based on the measurement of the longitudinal resistivity obeying 
the relation p^x = pix^~^^^'^^ 
increase of B while keeping 

V ~1 | |47[ |. If the number of reversed spins K increases, the effective Zeeman contribution jg/islBS is influenced via 
S, and the derivative of this expression with respect to B allows the determination of S. Schmeller et al. found 7 
reversed spins at g ~ 0.01, which is close to the theoretical prediction. Another modification of this technique can 
be realized by exerting hydrostatic pressure, which allows to change continuously g down to zero and to observe the 
increasing influence of spin-flips on the activation gap |138|. Despite these convincing expe rime ntal results, there is 
still some controversy about the influence of LL-mixing on the stability of these excitations |16C|. 

So far, we reviewed the theory of the low- lying charged excitations at = 1. Each of these two elementary 
excitations forms by definition also the ground state of a system with an area that contains one more or one less fiux 
quantum than at filling factor one. Since for small, but finite deviations from v — 1 the density of skyrmions is small, 
one can assume that these particles arrange for energetic reasons in a crystal similar to the formation of a Wigner 
crystal in the low filling factor limit of a QH-system. There are a few predictions on this Skyrme crystal. Except 
for very small deviations from \v ~ 1| ~ 0.04, where a triangular ferromagnetic arrangement of skyrmion centers is 
lower in energy, a HF-calculation favors a square-lattice arrangement of antiferromagnetically oriented spins at the 
skyrmion centers |161, 162 1. Any of the Skyrme crystal ground states is characterized by the long-range order in 
the spatial arrangement of the skyrmion centers as well as by the azimuthal orientation of the in-plane component 
of the local spin magnetization near each skyrmion center. Good agreement of a sp in magnetization calculation in 
HF-approximation with the results fro m NMR-experiments by Barrett et al. ]152| supports such a picture up to 
deviations ji' — 1| — 0.2, c/. Fi g. |20| |l6l|| . Moreover, a peak in the low-temperature specific heat can be related to the 
melting of the Skyrme lattice [163, 164 1. If the skyrmion density further increases, a transition to a fluid state seems 
to be reasonable. Recently, the ground state energy of a spin-polarized and an unpolarized liquid was calculated in 
the filling factor range 2/3 < < 1 [165|. It turned out that below v ~ 0.92 the spin- unpolarized liquid is favored over 
the square-lattice Skryme crystal, which questions the existence of a Skyrme crystal in such a wide range of filling 
factors. In any case, a very rich T — v — A^ ph ase diagram is expected to occur, which comprises various quantum 



and thermodynamic phase transitions [162 
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C. Exact neutral one-spin flip excitations a.t v — 1 



Thus far, we know the elementary charged excitations at fractional filling factors for spin-polarized systems and 
those for a model with spin at v — 1. These charged excitations become the constituents of neutral excitations 
above the ground states. They are either a well separated neutral quasihole-quasielectron pair or a neutral skyrmion- 
antiskyrmion pair far apart. In the Landau gauge, the quasihole-quasielectron pairs are the k oo limit of a whole 
class of excitations that describe one-spin-flip excitations above the ground state of a filled majority spin band 

at 1/ = 1. In a finite system with N = Ng, particles, N'^ states 

^E^""'"'-L^,+..tI*-i)' (93) 
g 

can be constructed. The summation over q = 2TTl/Ly {I = 0,±1,...) comprises all momenta in y-direction. The 
operator acting on IvPi/^i) is proportional to the spin lowering operator 5~(k), where the bar denotes the projection 
onto the LLL, and ^"(k) = J d^r e^'^'^^S^ {r) — Sx(ii) ~ iSyCk) is the Fourier transformation of the projected spin 
density 5"(r) = *j(r)^'^(r). Then 

5-(k) = e-^^^*'^^'^"^^ ^e--'^^^c et^e,^^^^ , (94) 

g 

and (|9^ ) can be written except for an irrelevant phase factor as 

|k) = ^^-(k)|*..i) . (95) 



The form of ( p3| ) ensures the normalization of the states |k). One can check that these states |k) are exact energy 
eigenfunctions of the LLL Hamiltonian in Landau gauge, cf. Eq. (p^, 

H = -\A.iN,-N,) + ^ Y: Wiq,p-p)cl^cl,^^,c^,^^^^,c^_^^^ . (96) 

PrP ,g 

Thus, the excitation energies esw(M) above the spin-polarized ground state |\E'i/=i) are described in the thermodynamic 
limit by the continuum of excitations 

esH'(k) = (k|i/|k) - 

= A, +A(a(0)-a(k)) . (97) 



The quantity 



a(k) = / ^y(q)e-^e'<i-(-x'')^l (98) 



is already known as the band formed by the eigenvalues of two attracting particles in the LLL, see Eq. (p8|). The 
three terms can be easily interpreted. First, is the one-particle energy that is necessary to promote an electron to 
the minority spin band. Second, Aa(0) is the exchange energy that is lost due to one electron less in the majority spin 
band. It is equivalent to the gross quasihole energy e_(l) in the majority spin band, see the Coulomb entry in Tab. ||. 
The third term, Aa(k), describes the attractive interaction between the electron in the minority spin band and the 
hole in the majority spin band and the energy equals that of two spinless, but differently charged fermions attracting 
each other, cf. Eq. (|2^). The two-particle complex of an electron and a hole can also be viewed as a spin-exciton or 
magneto- exciton, where the valence and conduction bands are identified as the dispersionless majority and minority 
spin bands of the LLL. 

Eq. ( p7|) describes for \h\£c ^ 1 the single particle character of the excitations, namely the interaction between 
the minority electron and majority hole whose distance is given by |k|£^. This can either be seen from ( p8| ) or by 
expanding ( ^7|) for large k£c- For example, in case of an unscreened Coulomb interaction, this yields 

eswi^) = A. + A Uo) -±\+o(^) . (99) 
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The second term is the energy of the particle-hole pair excitation Age-g/i at infinite distance and the lowest order 
correction due to the attractive Coulomb interaction between charges ±e of different sign at dimensionless distance 
kic- It is Aqe-qh — Aa(0) = X^/nj2. This behavior is analogous to that for the large distance limit of attractive 
quasiparticle interaction at = l/q, where the leading correction reflects the fractional charge of the quasiparticles 
by a prefactor v^. Below, in our diagrammatic analysis, we will see that it is necessary to work with a screened 
interaction a{k) instead of a{k). 

Of course, corrections to the spin-gap enhancement occur when one includes higher Landau level contributions. 
Then, the expansion parameter is \/{TiuJc) oc I/a/B, which vanishes in the LLL approximation. Besides analytical 
investigations, which are based on the leading order result of a perturbation theory ]135| ] and an RP A-ca lculation [167|, 
the reliability was checked by means of a variational quantum Monte Carlo (VQMC) approach ]16C|] . In that case, 
variational wavefunctions for the ground state as well as for the quasiparticles with a spin flip were constructed by 
multiplying these wavefunctions at = 1 by a Jastrow function Y[i<j e"^''^-* mixing in higher LLs. The perturbation 

theory by Sondhi et al. [135| leads in zeroth order of which can be written as the linear contribution in the 

dimensionless parameter Vs = ro/as, to a, particle-hole gap Aqe-qh with 



0.58 

7^' 



(100) 



In the high field limit 



0, we recover the LLL result, but for realistic values at = 1 and an electron density 



10^ 



it is in GaAs 



1.8. The Monte Carlo approach results in a reduction of the LLL gap of about 



45% to a value of approximately 0.7A compared with the LLL value 1.253A. Thus, in this parameter range the MC 
results are more reliable than the perturbation theory, which overestimates the influence of mixing for values larger 
than Ts ~ 1.2 leading to a value Aqe-qh — .515A. In summary, the VQMC interpolates between the perturbation 
results reliable at small and the RPA [167 , which is particularly suited to large values of the electron parameter 



The experimental determination o f the spin gap at = 1 by tran sport measurements of the activation gap [168|, 
magneto-capacitance measurements [169 , and optical methods [170[ do not support this ideal picture: first, the gap 
experimentally found, is much smaller, and, secondly, a linear instead of a square root dependence on B is measured. 
The main effect is attributed to the influence of disorder, but theoretical investigations could not support such a 
conjecture so far. Note, however, that these studies did not aim at the detection of skyrmion-antiskyrmion pairs. 

For small values [k[^c ^ 1, *• e. in the long-wavelength limit, the collective character of the excitations are of 
spin-wave character, and an expansion of (pT^) up to k^ yields 



esM/(k) = A.+AnpJlk^ 



(101) 



limfc_oA(a(0)-a(k))/(4<p^ 



The spin-sti ffness o 

cf. Eq. (pa), [135, 171 [, has its microscopic origin in the electron-electron interaction. Note th at in the limit k 



which appears as a phenomenological constant in field-theories, 

the 



excitation energy is that of the non-interacting system in accordance with Kohn's theorem |17S[. There are various 
other ways to derive the dispersion relation (pT\). Among these is the thermodynamic Green's function formalism 
where one has to calculate the spin-spin correlation function in the LLL. The identiflcation of the GF's pole in the 
limit T ^ leads again to the dispersion relation of Eq. (97). 

The excitations of Eq. (^) describe for each k bosonic excitations above the spin-polarized ground state since 
S = Sz = N/2 — 1. In systems exhibiting saturated ferromagnetism, such low- lying excitations are known as magnons 
[174[. Therefore one is inclined to explain the thermodynamics of the QHF at = 1 on the basis of non- interacting 
bosons, see Subsect. HID. However, such an approximation becomes more questionable if the number of these 



excitations is too large and the interactions between the bosonic spin- waves has to be taken into account. Such an 
effect can be seen from the diagonalization of a finite number of electrons on a sphere at v — 1 identifying one and two 
spin-wave states [173, 19S| [. Consequently, spin-wave spin-wave interaction plays a role in thermodynamics, especially, 
when the temperature increases and therefore limits the validity of a free magnon theory. 



D. Elementary theories of spin magnetization around v = 1 



The knowledge of the ground state of a QH many-particle system allows the calculation of the spin magnetization 
M{i',T = 0)) in dependence on the filling factor. But the investigation of the thermodynamic behavior of M{v,T) 
over the whole temperature range makes information about all excited states necessary. In general, such a task is too 
ambitious to be rigorously solved, and one needs approximative schemes to calculate thermodynamic quantities in a 
certain temperature range. 
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FIG. 18: Com parison of the exact dispersion curve (solid line) of Eq. @ with the spin-wave approximation (long dashed 
curve) in (101). Note the small shift at A: = due to the non-zero Zeeman term Az — 0.016A in comparison with the spin- wave 
bandwidth A5(0). A static screening wavevector k^c ~ 0.01^^^ is used. Around ki^ ~ 1, the two curves bifurcate indicating 
the diminishing collective character of the one-spin flip excitations fo r k£ r > 1 and the emergence of the minority spin-down 
electron-majority spin-up hole pair of single particle character. From [|l73|]. 



Due to the two-dimensionality, spontaneous magnetization at vanishing Zeeman term can only occur in the ground 
state, and we saw that the existence of a QHFM at = 1 is caused by interaction. On the other hand, a non-zero 
Zeeman term destroys the SU (2) symmetry and allows also for finite temperatures a non- vanishing magnetization. 
Before we study the influence of the interaction on the temperature dependent spin magnetization, we set out from 
the temperature and filling factor dependent spin magnetization of non- interacting electrons, which is given by 



M^'^Xi^, T = 0) = Mo (^6(1 + - 1)(^ - 1)^ 



(102) 



Eq. ( 102 ) shows constant magnetization as long as the majority spin band is filled up, but starts to decrease with the 
occupation of the minority band and vanishes for the completely filled lowest orbital LL at = 2. The generalization 
to arbitrary temperature yields 

^ ' 2 z + cosh(/3A J2) ' ^ ' 

where the fugacity z = e^^ (/i - chemical potential) is related to the filling factor v by 



z{v, T, A,) = ( \j (1 - cosh^ ^ + ^.(2 - i.) - (1 - cosh ^ ) , (104) 

see Fig. |l9|. For filling factor v — 1, we get a closed expression due to the independence of the chemical potential of 
temperature and Zeeman energy so that /i = and 

M(°) (i/ 1, T) = Afo tanh(/3A^/4) . (105) 

Note that the magnetization in a QH-system drops stronger than for a localized spin 1/2 system corresponding to a 
Heisenberg model on a lattice with zero coupling (J = 0), where M/AIq = ta,nh{f3Az/2). This is a consequence of 
the existence of spin as well as charge degrees of freedom, i. e. of the itinerant character of the QH-system. Indeed, 
an electron excited into the minority spin band can concomitantly reverse its spin and occupy any of the — N 
orbital one-particle states. In total, in a finite system iV^ excitations above the ground state are possible. 

From our previous discussion, we know that the zero temperature magnetization at = 1 is exact since all spins 



are polarized and show therefore saturated ferromagnetism. However, the non-interacting result of (102) is in obvious 
contradiction to our findings on the influence of the charged spin texture excitations, the skyrmions, on the ground 
state magnetization vs. shown in Fig. 1^. From a phenomenological point of view, the creation of one neutral 



quasihole at filling factor one changing t 



le filling factor infinitesimally to 1 ~ Sv leads to a rearrangement of the 



ground state, which leaves in average A holes in the majority spin band, i. e. K = A — 1 spins are fiipped in average. 
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FIG. 19: The temperature dependence of the spin magnetization of non-interacting electrons in the lowest Landau le vel a t 
filling facto rs 1/3, 1, and 5/3 for — 0.016A. The inset shows the ground state magnetization as a function of u, cf. Eq. (102). 
From |p|]. 



Analogously, the creation of a neutral quasielectron leads to K = S reversed spins in the minority spin band. The 
values A and S depend on the effective gyromagnetic factor g = A^/A. Thus, the spin magnetization about h' — 1 
can be expressed as 



M(i/,T = 0) = Mo 



9 9S ' 

6(1 - -A)-il- 2A)) + e{iy - 1)(— + 1-25) 



(106) 



This relation breaks down when the underlying assumption that independent skyrmions and antiskyrmions, respec- 
tively, form the ground state becomes invalid. Eq. (|102|), which is actually only true if 3> A, i. e. when the 



one-particle character dominates the magnetization curve, can be recovered when setting ^ = 5 = 1 in (106). For 
any value A and S larger than one, the magnetization drops precipitously from the maximum polarization at v = I. 
The deviation of the relative magnetization AI/Mq from the maximum value one is therefore up to second order in 
small values of 1|: 2{A- l)(|i^- 1| -I- 1| +■ . .) for < 1 and 2S{\iy- 1| - Ip -|- . . .) for > 1, respectively. 
Recent NMR-experiments by Barrett et al. [ |152| measure the Knight shift of '^^Ga, which is due to the additional 
magnetic field caused by the spins of th e 2DE S. This shift corresponds to the spin-polarization per particle around 
V — 1 at low temperatures, see Subsect. IIIF. Comparis on w ith the value of the left and right first derivative of the 
experimentally found spin magnetization curve at ^/ = 1 (|106| ) allows the determination of A and 5, respectively, see 
Fig. ^ Moreover, even the stronger drop in Barrett's curve for v < 1 breaking the symmetry of the magnetization 
about V = 1 in leading order of \y — 1| can be explained within this simple theory. The determined value from the 
curve are .4 = iS = 3.6 ±0.3 for an experimental ratio A^/A = 0.016. These results corroborate strongly the existence 
of finite-size skyrmions in a QH-system near v = 1. The found values are in fair agreement with theoretical estimates 
that predicted an effective spin of about 3.5 for this value of g |153 . The experimental values can be fitted within 
this simple theory for 0.9 < v < 1.1 and give a hint of the range of validity of the independent skyrmion assumption. 

A more ambitious theory to explain the ground state magnetization near v = 1 is bas ed o n the assumption of 
a Skyrme lattice, which we mentioned in conjunction with the discussion of skyrmions |161|. Results from HF- 
calculations are shown for a square and a triangular lattice in Fig. In the following, we will focus on the spin 
magnetization at exactly v = 1 whose study is essentially facilitated by the knowledge of the spin-polarized ground 
state. This is the reason why more ambitious theories on the temperature-dependence of the magnetization for filling 
factors unequal one, except for some FQHE filling factors, are still missing. 

Before we start to develop a many-particle theory of the magnetization at = 1, let us write down the spin 
magnetization for the case that the only excitations are those non-interacting bosons, which we described as one-spin- 
flip excitations in the last Subsection. Starting from the expression for the temperature dependent spin magnetization 
whose saturated magnetization value at zero temperature Mq is diminished by the number of occupied bosonic spin- 
waves with dispersion esw{^) 



M{T) = Afo - 5Mi3 "B(esw(k)) 

k 



(107) 
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FIG. 20: The ground state magnetization near v = 1 for Skyrme lattice states at ty pical values of g = A^/X. The open (7.05 T) 
and closed {B = 9.39 T) circles are the experimental results of Barrett et al. [152]. The legend indicates two different Skyrme 
lattice states: the SLA state is a square la ttice state with opposing skyrmion orientations; the TLF state is a triangular lattice 
with aligned skyrmion orientation. From [161|. 



we get after approximating the dispersion relation by Eq. (101) 



M(r) = Mo 



1 



In 



(1 



(108) 



This result is formally equivalent to the expression one gets for a localized spin system described by a Heisenberg 
model with spin S in an ex tern al magnetic field of strength Az/lgnsl when ksT <^ 2ttSJ, i. e. for low temperatures 
and spin S = 1/2. Hence, ( |108| ) represents the two-dimensional (linear T-term) analog in an external magnetic field 
of the T^/^-Bloch law in three dimensions and suggests that for low temperatures the thermodynamics is governed by 
independent spin- waves |174|. The paramagnetic behavior for temperatures greater than zero in the limit of = 
is signaled by the divergence of t he l ogarithm in the expression. 
For two temperature regions, (|l08|) simplifies to 



M{T) = Mo 



M{T) - Mo 



1 _ .^g-/3A, 

A-nps 



Attps 



A, 



A^ < keT < 47rp^ 



(109) 



In the first case, the magnetization is governed by activated behavior for temperatures much smaller than the activation 
threshold A^ due to the non-zero Zeeman term, the second case describes an approximately linear regime. The curve 
for the independent magnon approximation is shown in Fig. |2^. At temperatures above 0.14A/fcB, the magnetization 
becomes negative and the approximation breaks down. There are some obvious deficiencies, which are well known 
from the study of localized spin systems in three dimensions without magnetic field |176| . First, the number of 
fermionic excitations is restricted by the Pauli principle, while the occupation number of a bosonic state |k) can 
exceed one. Second, the interaction of spin-waves, whose existence was shown above, becomes at higher temperatures 
more important as the number of excited spin-waves strongly increases. This leads inevitably to an overestimate of 
the decrease of the magnetization and to negative values of M. Despite the restricted value of this simple the ory, it 
will give us some guidance for the low-temperature behavior of our many-particle approach discussed in Subsect. [II E . 



E. Hartree-Fock approximation and beyond: electron scattering on spin-waves 



Before we outline our many-particle theory taking into account the spin-wave excitations, let us analyze the con- 
sequences of the self-consistent Hartree-Fock approximation (SHF) at finite temperatures. Obviously, the SHF of 
Eq. (|o|) reduces to an algebraic equation for to be solved numerically for a given filling factor u = v^j + 1/^. With 
the particle- hole symmetry at v = 1, the temperature dependent SHF energies = e^^ — p satisfy the relation 
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FIG. 21: The temperature dependence of the spin magnetization at = 1 in the self-consistent Hartree-Fock approximation 
(dashed line), see below, and the result assuming free magnons (long-dashed curve). It is Az = 0.0 16A. For comparison, the 
results for non-interacting electrons (solid line, see also Fig. ^|) and of Barrett's experiment [152| (filled circles) are added. 
From 1^]. 



(,^^(T) = —£.^^(T) < 0, see Fig. Note that the chemical potential in the SHF at = 1 remains independent 
of temperature and Zeeman energy and is given hy ~ — Aa(0)/2. The temperature dependent energ y di fference 



— ^f^^ is the exchange splitting that shows oscillating behavior in dependence on the filling factor |177]. Here, 



we are rather interested in its temperature dependence for fixed = 1. The gap starts at zero temperature from 
its maximum value 2^^^ = -I- Aa(0) that is clearly dominated by interaction and reaches for T cx) the bare 
gap value reflecting the non-interacting situation. Experimentalists relate this enhanced spin gap to an effective 



g-factor, which can be determined from activation measurements [168|. However, comparison shows that the values 



found theoretically are too large and need to be corrected by incorporating t he in fluence of the finite width of the 



sample, screening effects, Landau level mixing as well as disorder, cf. Subsect. HI F 
The spin magnetization in the SHF yields 

M"^{T, iy = l) = Mo{v^^ - vf^) = Ma tanh(/3ef ^/2) (110) 

and is shown in Fig. Not surprisingly for a mean-field theory, the SHF exhibits in the limit A^ — incorrectly 
spontaneous magnetization due to a continuously broken symmetry in contradiction to the Mermin- Wagner theorem 
prohibiting this in dimension two for temperatures larger than zero. For non-zero A^, this leads to a softening of 
the magnetization curve, and the inclination point is a remnant of the mean-field phase-transition. In any case, 
also for non-zero A^, the magnetization is much too large when compared with experiment, see Fig. and such 
a behavior will prevail as long as the interaction is the dominating energy scale. Since the temperature dependent 
is proportional to the order parameter, the spin magnetization, the mean-field critical exponent (3 = 1/2 as 
well as all the other mean-field values can be found: a — for the specific heat, 7=1 for the divergence of the 
spin-susceptibility and S = 3 for the vanishing of the symmetry breaking field at the mean-field critical temperature. 
Formally, the discussion is analogous to the Weiss theory of ferromagnetism, where one finds a similar equation for 
the order parameter and whose results are independent of the spatial dimension. At temperatures above the mean- 
field critical temperature = Aa(0)/(4fcB), the vanishing magnetization leads to paramagnetic behavior, while 
for T < T^^ a thermodynamically stable solution unequal zero signals ferromagnetism. Note that the role of the 
symmetry breaking field is played by the Zeeman energy A^ and not by the external magnetic field. 

We see that HF theories are useful in predicting the ground state and the spin-flip excitations from the flUcd majority 
spin band to the minority spin band a,t v — 1. However, they yield rather poor results for the thermodynamic behavior. 
In particular, simple theories for the temperature dependent spin magnetization a,i v = 1 are insufficient: neither the 
non-interacting picture with its too strong drop of the magnetization at temperatures larger than the bare Zeeman 
energy A^ of about 2 K nor the SHF picture provides a trustworthy result when compared with experiment. In 
the latter case the exchange enhanced one-particle gap that only slowly decreases with temperature grossly 

overestimates the magnetization at temperatures below T^^ . Therefore, the magnetization curves of these two 
theories represent only lower and upper bounds for any improved theory as well as for the experimental data. The 
latter is even true if we keep in mind that in experiments additional corrections have to be taken into account, e. g. the 
finite width in z-direction. The SHF theory applied to the quantum Hall system aX v = \ has much in common with 
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FIG. 22: The self-consistent Hartree-Fock eigenenergy ^ (solid line) as a function of temperature at v = 1; = —(^^ 
because of particle-hole symmetry at u — 1. The magnetization M — Mo(i'| — vi) within the SHF is depicted for — 0.016A 
(dashed curve) and — O.OA (long dashed curve), respectively. Note the finite magnetization at low T in the latter case 
incorrectly indicating the existence of an ordered phase for T below Tc = o(0)A/(4fc_B). The uniform static inverse susceptibility 
is plotted as a dotted line in units of {16^Tl'^X)/{g^lBf for T > as well as for T < T^. From pTJ. 



mean- field theories of itinerant magnetism like the Stoner theory in the Hubbard model |178|. 

The improvement we propose accounts for the f act t hat spin-wave excitations are important for the suppression 
of the magnetization with increasing temperature |173|. We will systematically incorporate this mechanism in our 
diagrammatic treatment going beyond the SHF, i. e. we will identify those diagrams in an approximation of the 
self-energy E(w) that are responsible for such a mechanism. This approach is similar to an approximation that was 
applied to a single-band Hubbard model with strongly ferromagnetic ground states at zero temperature by Hertz and 
Edwards | p79| . 

Technically, we calculate the correction Scr(w) to the SHF self-energy so that S(j(a;) = + ^^(w). We 

express the self-energy by the four-scattering vertex F^''^, (1, 2, 3, 4), which is calculated in a ladder approximation 
exactly describing the interaction between a single electron in the minority spin-down band and a single hole in the 
majority spin-up band and vice versa. Since the pair-propagator Xo-.o-' enters the four-scattering vertex as a product 

of two GFs, the four quantities, i. e. the GF Q^, the corrected self-energy So-, the four-scattering-vertex r^"*^, and the 



pair-propagator Xo-.o-' form a closed set of equations, cf. Eq. ( 112 ). Note that the diagrammatic treatment calculates 
corrections to the SHF-GF and not as usually to the bare GF. The essential step is the summation of the particle-hole 



> 



q-q\kqll,h' 



FIG. 23: The self-consistent integral equation for the scattering vertex ^, evaluated in the particle-hole ladder approximation 
where a and o' =[. The first order term, the interaction vertex, has to be subtracted since it leads to a tadpole diagram 
for the self-energy, which has to be omitted due to the neutralizing background charge. 

ladder by solving the Bethe-Salpeter integral equation for the scattering vertex shown in Fig. p3| 

+ A(7/2 - 54, Ag; ^w„) = \W{q + Ag/2 - 94, Ag) 

/OO 
dq' W{q ~ q', Aq)ri% {q' + Aq/2 - q^, Aq; luj^) , (111) 
-00 

where the momentum part of the scattering vertex solely depends on the transfer of momentum qi — qi = q+Aq/2 — q4 
of one particle {q = (gi + q3)/2, Aq = qi — q^) and the difference of momentum between one outgoing and the other 
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ingoing particle Ag. 

Denoting T^^\, (iujn) as the four scattering vertex integrated over the two-dimensional momentum dependence of k 
our approximation can be summarized in the following set of four coupled equations 



(2^)H(l-a(k)Ax.,.'(^^„)) 



1 



Xa,a'{it^n) = " ^ ^ Scr (it'n (« + ^n) ) ■ (112) 

it'?! 

The self-consistent solution of this equation does not have stable solutions at low temperatures as was shown by 



Haussmann in conjunction with the problem of tunneling between two layers [137, 180 . Our treatment is based 
on an approximation, where the propagator in the pair-propagator and the self-energy is evaluated with the SHF 
Green's function (|79|). This leads to the correct low- temperature behavior as we see below. For example, the resulting 
self-energy correction for spin up I]|(w„) is 

' Jo 2 (iniyn + esw[k) - ^"^ ) 

Here, the Bose-Einstein distribution function nB{£sw{k)) of the temperature-dependent spin- wave dispersion 

-eswik) = A, + A(z.f ^ - i'f ^)(5(0) - a(fc)) (114) 

occurs. The expression for the minority spin can be found from the relation I]|(w„) ~ — S|(— Wn), which holds for 
1/ = 1 due to the particle-hole symmetry. At zero temperature, the correctness of the SHF is recovered since the 
self-ene rgy expression equals zero due to a vanishing numerator. 

Eq. ( |ll3[) and the corresponding minority spin expression are formally equivalent to the second order electronic 



self- energ y expression due to virtual phonon-exchange in a model of free band electrons coupled to phonons [181|. 
Eq. ( |ll3| ) describes the scattering of a majority spin electron into a minority spin state upon absorption of a spin- wave 
playing the role of the phonons. However, there is an important difference: the bosonic spin-waves carry spin S'z = — 1 
and not spin zero like the phonons. Therefore, it is not possible for the majority spin electron to have simultaneously 
scattering in the minority spin state Sz = —1/2 and emission of a spin- wave. Hence, an emission term as in the 
phonon-case is missing. Reversely, in the minority spin case scattering into majority spin states is only permitted 
when an emission of a spin-wave, but not of an absorption, occurs. In our case, the role of the electron-phonon 
coupling constant is played by the electron-spin-wave coupling, which is proportional to Xa(k). Note that, unlike 
the case of deformation potential phonon-coupling, the spin- wave coupling approaches for fc = a non-zero constant. 
This picture can be extended to the region of k£c 3> 1, where magneto-excitons with real space distance kl'^ can be 
absorbed and emitted, respectively. 

It is interesting to note that an unscreened Coulomb interaction leads to diverging self-energy expressions, which 
causes the slow asymptotic fall-off of a{k) ^ 1/fc, i. e. the long-range behavior at This can already be seen 

from the second order ladder diagram for the self-energy. Although screening does not play an important role in 
the physics of quantum Hall systems with gapped ground states at low temperatures, we are forced to introduce a 
screening mechanism. This is accomplished in our approach by replacing the unscreened interaction V{k) = 27r/fc by 
V{k) = 2TT / {k -\- ksc) , where we introduced the static screening wavevector k^c- Below, we will discuss the approximate 
temperature dependence of ksc accounting for the stronger screening with increasing temperature. Formally, any non- 
zero screening wavevector removes for real space distances ki'^ larger than 1 / kgc the divergence of the self-energy due 
to a quicker decrease of d{k). 

In summary, our partial summation of diagrammatic electron-electron contributions leads effectively to a model of 
SHF-electrons and spin-waves, which are coupled via an electron-spin-wave constant Aa(fc). It i s obvi ous that such a 



model goes beyond the naive picture of free magnons without coupling as described in Subsect. HID. 

The analytical structure of the retarded GF can be studied after analytical continuation to the real frequency axis 

G7\E) = i . . (115) 
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We are particularly interested in the properties of the spectral function A^(E), which is diagonal in the spin index 

A^E) = -2ImGT{E) (116) 

with the normalization condition 

dE , , , , , 

—A4E) = 1 . (117) 

The particle-hole symmetry at = 1 results in the relation A-^ [E] = A^ {—E) for the spin-dependent spectral functions. 

From the analytical structure of 'S^^*{E) we can infer that in the temperature dependent interval {£,f^, £,f^ — A^), 
the imaginary part of the self-energy is non-zero, and, hence, the retarded GF exhibits a branch cut in this energy 
range leading to a continuous contribution to the spectral function. This is due to the incoherent band contribution 
whose bandwidth {^f^ — Az) — — {I'f^ — i'^^)Xa{Q) agrees with the width of the spin-wave band. Outside this 
region, the zeros of the equation 

E-^^^ ^Ret^^^E) (118) 

define the quasiparticle poles in the spectral function. The self-energy has for T > always two quasiparticles poles, 
E^ and E'^ outside the spin- wave excitation interval. The former one is located below the spin-up SHF pole £,^^ ^ 
and the latter one in the positive energy region satisfying the relation E^ > — ■ 

This behavior is summarized in the explicit expression for the spectral function with the two quasiparticle terms 
and the incoherent band contribution 



A^E) 



2tt5{E-E-) 2TrS{E~E^ 



II 



dE \E=E-' 1-^ OE \E=E^ 



+ 1 



+ OiE-CWiC-E) i-^JruK-jE)) 

" {{E - - i?eE-*(i;))2 + (/mS-*(^))2) 

In order to get a qualitative feel for the low-temperature region, we discuss the behavior of the residues of the lower 
and upper quasiparticle poles for spin up 

= . (120) 

1-^ dE 1 



Dividing the self-energy integral (113) for S|(_E) into a collective contribution for k£c < 1 and a single-particle 



contribution when k£c > 1 shows the importance of the collective contribution for small k at low temperatures. The 



weight and location of the low-energy pole found by means of (|118|), (120) are 



= [l+£l / dkknB{esw{k))] ^ ~ 1 - / dkknB{esw{k)) 
Jo Jo 

= Cf^ - Aa(0)[l - z^"] , (121) 

where we require Airps S> ksT in order to keep the weight factor positive. It is remarkable that the coupling constant 
a{k) drops out of the expression for . The weight of the low-energy pole, which is one at zero temperature, is 
diminished by the number of excited spin-waves and its shift toward smaller values is proportional to the decrease 
of the weight. In the first line, we have assumed that the occupation number of bosonic spin-waves is small, an 
approximation that is in accordance with the omission of spin-wave interaction in our diagrammatic approach. A 
similar treatment for the upper spin-up pole yields 



1 

z+ = —-l = £l dkknB{£sw{k)) 
Zt Jo 



E+ = ^ - + \h{{))z+ . (122) 

This consideration shows how the high-energy pole is shifted from the upper limit of the incoherent contribution 
to higher energies. Thus, there exists repulsion of the two quasiparticle poles with increasing temperature. In the 
temperature range Anps ^ ksT ^ A^, it follows that 

+ ksT (kBT\ 
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while we get in the activated temperature region ksT <^ 



z+ = -^e "B^ . (124) 



These dependences lead to the same temperature behavior of the magnetization as described in (109). 

The discussion shows that in the low-temperature region the weights of the two quasiparticle poles exhaust almost 
the entire spectral weight. This remains even true at higher temperatures, where a transfer of weight from the lower 
to the higher pole occurs. This behavior is depicted in Fig. with results from numerical calculations of the spin-up 
spectral function for three different temperatures. For temperatures larger than ksT ~ 0.35A, the weight factors tend 
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FIG. 24: The spectral function A]{uj) at = 0.016A and screening wavevector kscic = 0.01 for temperatures ksT /\ — 0.05, 0.1 
and 0.2. The numbers indicate th e fra ction of the total spectral weight from the two poles and from the branch cut. The energy 
uj is given in units of X/h. From |l82|. 



to 0.5 from above and below, respectively |173|. The same happens for the Fermi factors at much higher temperatures 
ksT ^ Aa(0) so that the magnetization vanishes for T ^ oo as it should be. Although the incoherent spectral 
function contribution has only a ti ny in fluence on the spin magnetization, it becomes important when studying the 
nuclear spin-lattice relaxation rate |173|| . 



Before discussing the magnetization, let us comment on the role of screening [198|. As already mentioned, it is 
necessary to introduce an ad hoc screening mechanism. This adds to the bare Coulomb interaction Vc{k) a static, 
but temperature dependent screening wavevector ksc entering the Fourier transformation of the interaction as = 
27r/ {k-\-ksc)- At zero temperature, screening does not occur at any integer v like v = 1 due to the one-particle gap, but 
its importance grows with increasing temperature larger than zero. A dynamical screening t heor y in random-phase 



approximation down to intermediate magnetic field strengths was developed by Smith et al. | 167 |. Here, we restrict 
ourselves to the LLL and start from a Thomas-Fermi theory expression B 

The existence of incompressible states at integer filling and zero temperature means k ~ Q and ensures the vanishing 



of ksc- Eq. (125) can be written as 



ksc -^f^ ^^MEW{E){1 - nF{E)) , (126) 

which allows to use the calculated spin-wave spectral function. This is an implicit equation for kgc since the spectral 
function depends on ksc- Explicit determination of ksc was carried out by evaluation of the spectral function entering 
the r. h. s. of ( |126D , which depends on ksc via T^(k). We found for the two temperatures ksT = 0.09A and fc^T = 0.18A 
wavevector values ksc^c = 0.01 and ksc^c = 0.1, respectively, which shows that ksc is more than a simple fit parameter. 
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The knowledge of the spectral functions Aa- (E) , which is the many-particle analog of the one-particle density of 
states, allows the calculation of thermodynamic quantities. The temperature-dependent spin magnetization Eq. ( |7l| ) 
can be written for = 1 as 

M{T) = -Mo r ^tanh{^)A^iE) . (127) 

The one-particle results aX v — \ for non-interacting electrons and within the SHF can be easily derived. The usage 
of the spectral functions Ai^{E) = 2tt5{E — ^z/"^) and Ai^{E) — 2'k5{E — ^f^), respectively, yields the expressions 



for the magnetization of Eqs. (105) and (110), see Fig. kII. The temperat ure d ependence of the magnetization over 



t he e ntire temperature range can only be found by numerically evaluating ( |127] ) on the basis of the spectral function 



(119) 



F. Comparison of theory and experiment: the spin magnetization at v = 1 

In this Subsection we describe shortly two other theoretical approaches that also aim at the understanding of 
the physics at filling factor v = \^ before we introduce experimental techniques allowing the study of the spin 
magnetization. Eventually, we will extend our theoretical results by taking into account deviations from our model 
of an ideal 2DES and compare these results with those from experiment. 

To date, there exist two strategies to attack the problem of calculating the magnetization at = 1: first, by starting 
from the electronic many-particle system emphasizing the charge and spin degrees of freedom, and second, by mapping 
the original problem onto a ferromagnetic lattice Heisenberg model and deriving an effective continuum quantum field 
theory (CQFT). Examples of the former approach are a work by H aussmann p8U and our calculation discussed at 



length before. The latter was put forward by Read and Sachdev |184| on a mean-field level and later extended by Timm 



et al. [185| calculati ng co rrections. Numerical studies like exact diagonalizations in various geometries [144, 173, 186 



or MC calculations [187| play rather an intermediate role. Although they struggle in general with finite-size effects, 
they serve as a check of analytical theories. 

The theory developed by Haussmann is particularly successful in predicting the spin magnetization in the high- 
temperature range down to temperatures, at which the relative spin magnetization does not exceed a value of about 



0.4 [180 . It is a so-called modified self- c onsi stent random-phase approximation and was originally developed in order 



to explain tunneling experiments [136, 137 [. The idea is quite general and has its origin in the treatment of the 



independent boson model, an exactly solvable model that describes the interaction between one electron of fixed 



energy and independent bosons linearly coupled to the single electron [ 181 [. This set of ideas can be transferred to 
the problem of interacting electrons in the LLL using a diagrammatic formulation. Breaking the Coulomb interaction 
in two parts, where one part can be viewed as bosonic, allows to transform the original interacting problem into an 
interacting electronic problem with a residual interaction. Applying a modified self-consistent RPA he was able to 
determine the self-energy of the transformed problem, but also the spectral fmiction A„{E) of the original problem, 
which consists of two Gau ssian peaks each of them weighted by Fermi factors. This reminds of the similar structure of 
the spectral function ( |lig| ), where, however, delta peaks occur. Although Haussmann's theory becomes instable when 
calculating ^o- at low temperatures, the form of the spectral function seems to be rather independent of temperature. 
At zero temperature neglecting the Fermi factors, two delta-peaks of a distance that agrees with the exchange energy 
in the LLL were found. Due to the temperature suppression of the upper pole, only one pole as in the SHF occurs. 

Unfortunately, this theory pretends also a phase transition when ^ at fc^Tc — 0.041A. This is just the 
temperature region, where this theory breaks down. We think that this fact is due to the neglect of collective 
excitations as in the SHF-approximation. 

The field theoretical description is based on the observation that the quantum Hall ferromagnet at = 1 and a 
ferromagnetic quantum- mechanical Heisenberg model for spin S = 1/2 on a square lattice in an external magnetic 
field are phenomenologically equivalent because of a ferromagnetic ground state and the existence of spin-waves as 
low-lying excitations. The Heisenberg Hamiltonian is given by 

i? = -J^ SiSj-A,^5,,i , (128) 

<ij> i 

where < ij > means summation over nearest neighbor sites of a quadratic lattice. The one-spin flip excitations 

[k>=^ye*''''^57[0> (129) 

3 
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N - number of lattice sites, describe spin-waves of energy — —J{cos{kxa) + cos{kya)) above the spin-polarized 
ground state |0). The mapping to the original problem is accomplished by matching the coupling constant of the 
Heisenberg model J = Aps > to the spin-stiffness found in the long-wavelength limit ( |lOl| ) of the dispersion relation 
of the microscopic Hamiltonian, and which is proportional to the exchange energy Aa(0). For an unscreened Coulomb 
interaction, we find J/A = a(0)/(47r) = l/(4-\/27r) ~ 0.0997, while screening and finite width effects result in a 
reduction of the spin-stiffness. It is obvious that such a mapping neglects the charge degree of freedom and reduces 
the problem to a localized spin problem, whose quantum number k is restricted to the first Brillouin zone. Moreover, 
the neglect of itineracy leads to differences in the physically less interesting high-temperature region. The leading 
order in a high-temperat ure ex pans ion of the Heisenberg model is proportional to /3A^/2 and is twice of the value 
for the QHFM, see Eq. (105) |183|. Another basic question concerning the completeness of the model was raised 



when the infiuence of the Hartree term (86) on the magnetization was studied |188{ . As we have already seen, a 
systematic gradient expansion of the order parameter field ni(r) contains additional terms beyond the leading non- 
zero Heisenberg model contribution. Brey showed that the next order contribution, the Hartree term, does have a 
non-negligible influence on the magnetization in the temperature range 0.05 < ksT/X < 0.15 leading to an increase of 
the magnetization by approximately thirty per cent. This is due to the preference of a more homogeneous topological 
density in the Hartree term suppressing spin-fiuctuations. This conclusion was drawn from classical MC results, 
which should agree quite well in this temperature range with the quantum- mechanical result. Hence some doubt is 



left whether the Heisenberg model (128) can really serve as an a ppropriate model in such a wide temperature range 



To get immediately numerical results, one can evaluate ( |l2q ) by means of quantum Monte Carlo calculations on 



a lattice as it was done by Henelius et al. obtaining very accurate results [187|. An analytical treatment of the 
Heisenberg model is possible by generalizing the original model with SU(2) = 0(3) symmetry in spin space to the 
cases of SU{N) and 0{N) symmetry in a Schwinger boson representation. That means, when treating the general 
case of spin S, i. e. = S{S + 1), that the entries of the N x N SU{N)-ma.trix = bl^bp are composed out of N 
bosonic operators (Schwinger bosons). Additionally, they have to satisfy the constraint ^a^a = NS. Performing 
the continuum approximation leads to a CQFT whose partition function is evaluated in the N oo limit, which 
is known as the large N mean- field theory. Improvements are expected from the calculation of 1/A'^-corrections for 
N = 2 and 3, respectively, which are not really small. In total, this approach provides four different results, depending 
on the symmetry and on the accuracy of the large N calculation. 

The most accessible observable reflecting the peculiar features of this correlated electron system at = 1 is the spin 
magnetization. So far, three methods were applied for a more or less direct measurement of this observable at v — I. 
These are optically pumped nuclear magnetic resonance (OPNMR) measurements, the magneto-optical absorption 
technique and photo-luminescence experiments. 

Information about the spin magnetization of a 2DES in an electron doped GaAs-Al^Gai-^xAs quantum well can be 
extracted from NMR experiments. Essentially, a nucleus with spin I and magnetic moment /xat is not only Zeeman 
coupled to the external magnetic field B — Be^, but also to the electron's spin S of the 2DES via the hyperfine 
interaction, i. e. 

Hn ^ -gNf^NBI + AIS, (130) 

where A is the hyperfine coupling constant. Rewriting the coupling between nuclear and electron spin-operators as 
IS = (7+5^ + /^5+)/2 + IzSz allows to distinguish between the influence of static and dynamical mechanism on 
the NMR-spectra. The expectation value of the transversal part of the hyperflne interaction vanishes for an external 
fleld in z-direction. The additional magnetic field at the location of the nucleus, which is proportional to (Sz), causes 
a shift of the resonance NMR-signal with respect to the resonance in a GaAs-sample without hyperflne coupling. 
This shift of the NMR-resonance, the Knight shift, can be used to extract the spin magnetization of the 2DES from 



the measured NMR-spectra [189 . Unfortunately, the shift does not allow the determination of the absolute value 
for the spin-polarization, however, there exist often reference points. For example, the maximum Knight shift at the 
lowest available temperature can be identifled with the occurrence of saturated ferromagnetism at exactly v = 1. 
Deviations from this ideal behavior appear as measurements at zero temperature are not possible, and the influence 



of disorder diminishing the spin-polarization is often unknown [190|. Hence, some source of uncertainty remains when 



evaluating experimental data. Such an NMR-technique was developed and applied to the 2DES in a strong magnetic 



fleld by Barrett and collaborators [g6|, 152]. One of the main obstacles was to reach a sufficiently large NMR-signal 
from the ^^Ga-nuclei. An NMR-signal enhancement of ^ 100 could be obtained by means of an optically pumped 
NMR-technique (OPNMR). Moreover the technique allows to distinguish the signal from the nuclei within the well 
from that of the nuclei in the barrier region p6[ | . More information can be obtained from the time-dependence of the 
Knight-shift, which is related to the spin-lattice relaxation rate l/Ti of the out-of-equilibrium nuclear spins due to the 
coupling via the transversal terms to the electrons. The measurements were done with GaAs-Alo,iGao,QAs multiple 
quantum wells of width 300 A separated by 1800 A wide barriers exhibiting typical electron densities of around 
lO^^cm"^ for the 2DES in the well and a high mobility larger than 10^ cm^ /Vs. The strength of the magnetic flelds 
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was between 3 and 12 T in dependence on the investigated filling factor. I nitia lly, the lowest available temperatures 



were about 1.55 K, but meanwhile values smaller than 0.3 K are possible [159| 



The magneto-optical method initiated by Goldberg's group 191 1 utilizes the temperature dependence of the 
absorption coefficient a{ijj,T) upon shining light onto a Ga^s-heterostructure. The peaks in the absorption spectrum 
can be associated with optical transitions of electrons from an initial state i in one of the valence bands into a final 
state i in the LLL with spin a. Integrating over the entire peak w. r. t. the frequency of the experimentally determined 
absorption curve assuming some frequency cut-off yields the intensity lij = J dway (w), which is proportional to the 
optical matrix element (uj) for this interband transition as well as to the total number Na^ = — N„ of available 
or non-occupied states in the LLL with spin a, i. e. 

h, = CUjNa^ . (131) 

A theoretical calculation of the fij and the determination of the factor C, which is approximately a constant indepen- 
dent of the transition and temperature, yield the quantities Na^ and Na^ ■ These are related to the spin magnetization 

by 

Mm^M.^^M,^^^^^ (132) 

The unknown total particle number N in the LLL can be inferred from the relation {Na^ + Na^)/N ~ (2/j^) — 1. The 
advantage of this method is the determination of the absolute value of the spin magnetization, while, on the other 
hand, additional theoretical calculations of the optical matrix elements are necessary. 

Theoretical simplifications made so far comprise mainly the omission of the finite width w and the finite barrier 
height of the quantum well containing the 2DES, the neglect of higher orbital Landau levels as well as of disorder. 

The first issue can be taken into account by a form factor i^(k, w) leading to an effective interaction I4//(k, w) = 
F (k, w)V (k) , which is due to the extension of the envelope wave function (f>o{z) in z-direction instead of a delta- 
function for zero width. Assuming an infinite barrier w ell h eight and a symmetric charge distribution p(z) = |(;/)o(z)p — 



9{w — z)(2/w)sin {ttz/w) as in Barrett's experiment |152|, we obtain for the form factor 



327r4(e''="' - 1) Btt^ 3kw 

^ ''^>~ (fcw(47r2 + k^w^)y ^ (fcw(47r2 + k^w^)) ~^ {An^ + k^w^) ^ ' 

altering the short-range behavior for kw 3> 1 . In the phenomenological field-theory, the spin-stiffness ps becomes 

= dkk^Veff{k,w)expi-k^£l/2). (134) 
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For a typical width of 30 nm and fields of i? = 7.05 T as in Barrett's experiment, the width in units of the cyclotron 



length is w ~ 3.ll£c- The spin-stiffness becomes then approximately half of its zero-width value |173| and causes a 
considerable decrease of the energy scale. 

Since the ratio of Landau level distance and interaction energy scales with the square root of B, sec (pT[), the 
correction due to higher Landau levels depends on the magnitude of the magnetic field used in experiment. In 
Barrett's experiment, there is X/{fiuJc) = 0.966, which is not really much smaller than one as it should be. In fact, 
allowing arbitrary values for the magnetic field, but still keeping v = 1, the task becomes much more complicated as 
even saturated ferromagnetism cannot be simply assumed, cf. an example at ly — 2 with a first-order phase tr ansi tion 



from paramagnetism to ferromagnetism upon lowering the magnetic field |192]. The Monte Carlo results |16C|| for 
the one-spin flip gap indicate a considerable decrease in the k — > oo-limit of the spin- wave dispersion. Nevertheless, 
it remains difficult to estimate the quantitative influence of a renormalized dispersion relation on the temperature 
dependent magnetization. In any case, the magnetization should be diminished due to the consideration of LL mixing. 

The largest uncertainty arises from the influence of disorder. Almost all theories aiming at the calculation of 
thermodynamic properties neglect it. In spite of the high-mobility of the samples, disorder occurs, which leads in 
general to a broadening of the Landau levels. This can be taken into account in a crude approximation by the 
self-consistent Born approxi matio n (SCBA) resulting in the somewhat unphysical semi-elliptical broadening of the 
one-particle density of states |177[ . Similarly, t he m any-particle gap of the FQHE is reduced, what was quantitatively 
shown for short-range impurities at z/ — 1/3 |l93| , as well as the disappearance of the FQHE when decreasing the 
mobility of the samples below a certain value |194[ . Qualitatively, disorder should decrease the magnetization of the 
ground state due to reversed spins in regions of strong potential fluctuations. 

Comparison of theoretical and experimental results is for many reasons not easy. Simplifications and approximations 
in theory as well as non-unique experimental conditions influence the results. Nevertheless, the 2DES in the QH-regime 
offers outstanding conditions for comparison of theory with experiment in a highly-correlated system. Therefore, let 
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FIG. 25: Results for M{T) at A2 = 0.016A for various values of the screening vector ksc and the widths w and a, respectively. 
The quantity a used in finite-size numerical calculations describes the width of a Gaussian charge distribution in ^-direction 
instead of the width of a hard-wall quantum well used elsewhere. 1. free electrons: w = - dotted curve. 2. SHF: 
ksc = 0.01^7^, ui = - long-dashed line. 3. our theory: ksc = O.Ol^^^i^f ~ 3.11^c - solid line. 4. our theory: ksc = 0.1^7^, = 
S.ll^c - solid line with dots. 5. exact diagonalization on the sphere: A*' = 9, ksc = G,a = Q solid line with cross es. 6. exact 
diagonalization on the sphere: N = 9, ksc = Q,oi = 2ic - soli d lir ie with circles. 7. Barrett's experimental data [152|: w = S.ll^c 
- filled points. 8. 0(A'^)-field theory with l/A''-corrections |18!:]: w = S.ll^c - long dashed line with squares. 9. SU{N)-iie\A 
theory with l/A''-corrections ]l85| ]: w = S.ll^c - long dashed line with crosses. 10. Monte Carlo results for the Heisenberg 
model |185|: w — 3.11£c - dot-dashed line with stars. From 173|. 



us quickly compare various theoretic al an d experimental magnetization results, which are depicted in Figs. p5| and p^, 
for a more thorough discussion, see [173|. Our comparison with experiment is mainly based on Barrett's data, while 
results from magneto-absorption experiments are only shown in Fig. The SHF and the free particle results serve 
as upper and lower bound for any magnetization data. Barrett's NMR data are in between these two curves except 
for temperatures below 0.02, where they seem to be too high, and for ksT/X > 0.09, where they even fall below the 
free particle curve. A possible explanation for the former behavior is the measurement at ly — 0.98 instead of = 1 
as well as disorder preventing full polarization at zero temperature and therefore keeping the curve smoother when 
approaching zero temperature. Not surprisingly, Goldberg's data whose technique allows the determination of the 
absolute value of the relative spin-polarization are smaller in this temperature range, see Fig. ^ 

Comparison of the numerical data from the MC calculation and the exact diagonalization results on the sphere 
for nine particles assuming a finite width gives additional important information. The MC data are always larger 
than the diagonalization data. While this has to be true at larger temperatures due to the omission of the itinerant 
character of the electrons in the Heisenberg model, the discrepancy at lower temperatures is mainly due to the strong 
finite-size corrections of the diagonalization data at small A^. The temperature, at which these finite-size corrections 



become important, can be estimated from numerical results for the susceptibility on the sphere |173 . Therefore, the 



MC data seem to be more trustworthy at low temperatures. The results of the CQFT partially differ in dependence 
on the approximation used. At low temperatures, the SU{N) mean-field theory is most appropriate as the free 
magnon scenario (108) can be recovered. Except for tiny deviations, this is also true for the SU{N)-theoiy with 
1/A^-corrections. At temperatures above 0.4, the curves start to diverge. The mean-field curve is above the MC 
data while the 1/iV-correction curve drops so strongly that it becomes negative similar to the result of a free magnon 
theory of Fig. B^. On the other hand, the 0{N) results do not yield the correct free magnon behavior, but agree 
well with the MC data. Here, the 1/iV-corrections represent a notable improvement over the mean-field 0(A^)-theory. 
In spite of doubts regarding the completeness of the model, the results are impressive for the experimental value of 
g. These findings within the CQFT corroborate the picture that the collective behavior dominates in the considered 
temperature range for 47rps ^ A^. Such an interpretation gets further support when comparing the dispersion curves 
for the QHFM and the Heisenberg model, which start to diverge in a system of finite width only at one-spin flip 
excitation energies that correspond to temperatures of more than 50K, cf. Fig. 

Among the many-particle approaches, our theory is evidently an essential improvement over the SHF-theory. The 
correct low- and high-temperature behaviors are additional pleasing features. Moreover, the discussion of the low- 
temperature excitations showed that in the case of vanishing Zceman term the spin-wave excitations suppress any 
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FIG. 26: C omp arison of the temperature dependent spin magnetization M{T) from magneto-absorption data [191] with Bar- 
rett's data 1 195 1 . The results for non-interacting particles with Zeeman energy = 0.016A and our many-particle result are 
also shown. Note the higher magnetization at higher temperature from magneto-absorption experiments, which do not fall 
below the free particle result as Barrett's data do. 



spontaneous magnetization at temperatures larger zero. On the other hand, the deviations from experiment are 
still large at moderate temperatures despite improvement when accounting for the finite width of the quantum well. 



two 



The theory of Haussmann suggests that at higher temperatures screening plays an important role. In Fig. 
curves with different screening wavevectors k^c = O.Ol^J^ and are shown. They correspond to those two 

temperatures 0.09A/fcs a nd O .lSA/fcg, respectively, that we found from the self-consistent determination within the 
RPA-calculation, c/. Eq. ( |126D . Therefore, we expect the curve with the correct temperature dependence to interpolate 
between the two curves in the temperature range between 0.09 A/fcs and 0.18 A/fc^. 

Within our diagrammatic approach, skyrmions cannot be taken into account, but their influence on the spin 
magnetization is presumably not very strong. Although skyrmion-antiskyrmion neutral excitations are more effective 
in decreasing the magnetization due to the larger number of flipped spins, their thermodynamic weight is rather 
small because of their comparably large energies. Moreover, the experimental Zeeman term value A^/A = 0.016 lets 
dramatically shrink the energetic preference for these excitations over the neutral quasielectron-quasihole excitations. 

The knowledge of the spectral function allows also comparison with other experimental observables at filling factor 
one. For example, the filling fact or d ependence of the spin- lattice relaxation rate l/Ti, which was systematically 
studied at temperature T — 2.2K [195| shows a deep dip aX v = 1 due to the lack of low- lying excitations in contrast 
to those existing above the ground state of a Skyr me crysta l. Unfortunately, only single values for the temperature 
dependence of 1/Ti at e xactly v — 1 are available |195, 196|. On the basis of a Korringa theory of nuclear magnetic 



resonance in metals |18!;] employing an approximation without vertex corrections we predict a Gaussian like maximum 
i n the tem perature dependence of the relaxation rate |173|. The rate was also discussed in the framework of the CQFT 



[184, 187| 



A traditional method to measure spectral functions is tunneling |197| . In our case, a current can flow between two 
parallel 2DES layers in a magnetic field when a voltage V between the layers is applied. In the case of two weakly 
coupled layers, the current can be expressed as the convolution of the energy dependent spectral functions. We 
predict pronounced peaks in the /-^-characteristic if the voltage V equals the temperature dependent difference of 
the quasiparticle peaks, i. 



eV = E+ 



^1 - 



173 . Although tunneling experiments were performed in 



the range 0.2 < v < 0.9 for spin-polarized electrons |136|, experimental data at filling factor one are not available to 
date. 



IV. SUMMARY AND OUTLOOK 



This paper is concerned with the physics of interacting electrons in two-dimensions in a strong magnetic field 
neglecting the influence of disorder. Such systems can be realized in certain semiconductor structures of sufficiently 
high mobility, where the electrons form a 2DES. The electronic degrees of freedom can be well observed at low 
temperatures. There are some unique properties of such a 2DES. Among them are: the position of the chemical 
potential can be tuned, the band structure of a clean sample in a magnetic field is exactly known, and the coupling 
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of the electron's spin to an external magnetic field can be changed. These properties allow a much easier comparison 
of theory with experiment than in other strongly correlated systems. In the strong magnetic field limit, however, we 
encounter the common problem that any naive perturbation theory fails. This fact as well as the many unpredicted 
experimental observations in the quantum Hall regime have been a challenge to theory for more than twenty years. 

Our theoretical treatment is restricted to the physics of the lowest orbital Landau level based on a microscopic 
electronic Hamiltonian. First, we study the extreme case of an infinitely strong magnetic field, where the physics 
is solely determined by interacting spinless electrons. Although this case is rather artificial, it provides important 
information, which can explain the incompressibility of the many-particle ground state at certain filling factors. This 
is a necessary ingredient for the understanding of the fractional quantum Hall effect. In particular, we investigate 
the properties and quality of various quasiparticle trial wavefunctions. The inclusion of the spin degree of freedom is 
essential at and near filling factor one. We show that at i' = 1 the ground state exhibits saturated ferromagnetism, 
even in the absence of a symmetry breaking field, while novel spin-structures appear near = 1 if the effective 
gyromagnetic factor g is small enough. The thermodynamics of such a quantum Hall ferromagnet is studied at 
exactly filling factor v — I in the framework of a many-particle theory emphasizing the role of spin-wave excitations 
and screening. We investigate the temperature dependence of the spin magnetization in the presence of a small 
Zeeman term and compare with other theoretical approaches and with experimental data. 

Although predictions of future research are rather limited as many of the unprecedented experimental observations 
show, let us shortly speculate, what will be of particular theoretical interest in the next time. So far, we know 
quite a lot about the physics at single filling factors, where experimental effects are most pronounced. On the other 
hand, there remain large regions along the filling factor line, where much less is known. Incompressible, compressible 
metallic-like and insulating ground states are known at certain filling factors, but the transition between them is 
hardly understood. Support from the experimental side is rather restricted due to the occurrence of disorder. Thus, 
the consideration of disorder is an interesting and promising field. Unfortunately, quantitative control of disorder in 
experiment is limited. Moreover, such investigations could also diminish the gap between the rather loosely connected 
fields of interacting electrons in a magnetic field and non-interacting electrons in a random potential. 

Another interesting field seems to be the systematic study of higher Landau levels. Recent experiments suggest 
that some of the spectacular effects of the lowest orbital Landau level disappear when the electrons reside in higher 
Landau levels. Nevertheless, our knowledge about this region is rather sparse. In all, the study of the transition from 
the high-field region down to the low-field region is still almost unexplored. 

In closing, the two-dimensional electron system in a strong magnetic field exhibits a lot of different physical phases, 
and an astonishing number of theoretical concepts finds their application. In view of the progress of the last twenty 
years and the number of unsolved questions, future experiments and theories will certainly deepen and enlarge our 
current understanding. 

The present work is a shortened version of the Habilitation thesis of the author submitted to the Otto- von- Guericke-Universitat 
Magdeburg. Discussion and collaboration with W. Apel, A. H. MacDonald, J. J. Palacios, S. Barrett, S. M. Girvin, B. Goldberg, 
R. Haussmann, P. Henelius, C. Timm as well as the support by H. Bottger, P. Fulde, B. Kramer, and U. Rofiler are gratefully 
acknowledged. The work was partially supported by a grant of the DAAD. 
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